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Abstract High-energy X-rays and y-rays from solar flares were discovered just over fifty 
years ago. Since that time, the standard for the interpretation of spatially integrated flare 
X-ray spectra at energies above several tens of keV has been the collisional thick-target 
model. After the launch of the Reuven Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI) in early 2002, X-ray spectra and images have been of sufficient quality to allow a 
greater focus on the energetic electrons responsible for the X-ray emission, including their 
origin and their interactions with the flare plasma and magnetic field. The result has been 
new insights into the flaring process, as well as more quantitative models for both electron 
acceleration and propagation, and for the flare environment with which the electrons interact. 
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In this article we review our current understanding of electron acceleration, energy loss, and 
propagation in flares. Implications of these new results for the collisional thick-target model, 
for general flare models, and for future flare studies are discussed. 

Keywords Sun: flares • Sun: X-rays • Gamma rays • Sun: radio radiation 


1 Introduction 

A primary characteristic of solar flares is the acceleration of electrons to high, suprathermal 
energies. These electrons are observed directly in interplanetary space, and indirectly at 
the Sun through the X-ray, /-ray, and radio emissions they emit (Hudson and Ryan 1995). 
Understanding how these electrons are produced and how they evolve is fundamental to 
obtaining an understanding of energy release in flares. Therefore, one of the principal goals 
of solar flare research is to determine when, where, and how these electrons are accelerated 
to suprathermal energies, and what happens to them after they are accelerated to these high 
energies. 

A major challenge to obtaining an understanding of electron acceleration in flares is 
that the location where they are accelerated is not necessarily where they are most easily 
observed. The flare-accelerated electrons that escape the Sun are not directly observed until 
they reach the instruments in space capable of detecting them, usually located at the distance 
of the Earth from the Sun. The properties of these electrons are easily modified during their 
long journey from the flaring region to the detecting instruments (e.g., Agueda et al. 2009). 
Distinguishing flare-accelerated electrons from electrons accelerated in interplanetary shock 
waves is also difficult (Kahler 2007). 

The electrons that are observed at the Sun through their X-ray or /-ray emissions radiate 
most intensely where the density of the ambient plasma is highest (see Sect. 2). Therefore, 
the radiation from electrons in and near the acceleration region may not be intense enough to 
be observable. Although these radiating electrons are much closer to the acceleration region 
than those detected in interplanetary space, their properties can still be significantly modified 
as they propagate to the denser regions where they are observed. The radio emission from 
the accelerated electrons also depends on the plasma environment, especially the magnetic 
field strength for the gyrosynchrotron radiation observed from flares (Bastian et al. 1998). 
Therefore, determining when, where, and how the electrons were accelerated requires a 
substantial amount of plausible reasoning. 

Here we focus primarily on the X-ray emission from the accelerated electrons. Interplan- 
etary electrons and low-energy emissions are addressed by Fletcher et al. (201 1), while the 
/ -ray emission is addressed by Vilmer et al. (20 11), and the radio by White et al. (20 11). The 
X-rays are predominantly electron-ion bremsstrahlung (free-free radiation), emitted when 
the accelerated electrons scatter off ions in the ambient thermal plasma. Issues related to the 
emission mechanism and deducing the properties of the emitting electrons from the X-ray 
observations are primarily addressed by Kontar et al. (2011). Here we address the interpre- 
tation of the X-ray observations in terms of flare models, and consider the implications of 
the observations for the acceleration process, energy release in flares, and electron propaga- 
tion. Specific models for particle acceleration and energy release in flares are addressed by 
Zharkova et al. (201 1). 

The accelerated electrons interact with both ambient electrons and ions, but lose most of 
their energy through electron-electron Coulomb collisions. Consequently, the brightest X- 
ray sources are associated with high collisional energy losses. These losses in turn change 
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the energy distribution of the radiating electrons. When the accelerated electrons lose their 
suprathermal energy to the ambient plasma as they radiate, the source region is called a 
thick target. Electrons streaming downward into the higher densities in lower regions of the 
solar atmosphere, or trapped long enough in lower density regions, will emit thick-target 
X-rays. Hence, thick-target models are important to understanding the origin and evolution 
of accelerated electrons in flares. Thick-target X-ray emission is addressed in Sect. 2. 

The total energy carried by accelerated electrons is important to assessing acceleration 
models, especially considering that these electrons carry a significant fraction of the total 
energy released in flares. Also, the energy carried by electrons that escape the acceleration 
region is deposited elsewhere, primarily to heat the plasma in the thick-target source regions. 
The X-ray flux from flares falls off rapidly with increasing photon energy, indicating that the 
number of radiating electrons decreases rapidly with increasing electron energy. Therefore, 
the energy carried by the accelerated electrons is sensitive to the value of the low-energy 
cutoff to the electron distribution. The determination of this low-energy cutoff and the total 
energy in the accelerated electrons is addressed in Sect. 3. 

In the standard thick-target model, the target plasma is assumed to be fully ionized. If 
the target ionization is not uniform, so that the accelerated electrons stream down to cooler 
plasma that is partially ionized or un-ionized, the X-ray spectrum is modified. This is ad- 
dressed in Sect. 4. 

Observations of the radiation from hot flare plasma have shown this plasma to primarily 
be confined to magnetic loops or arcades of magnetic loops (cf. Aschwanden 2004). The 
observations also indicate that the heating of this plasma and particle acceleration initially 
occur in the corona above these hot loops (see Sect. 12.2 and Fletcher et al. 2011). When 
the density structure in these loops is typical of active region loops, or at least not highly en- 
hanced above those densities, the highest intensity, thick-target X-ray emission will be from 
the footpoints of the loops, as is most often observed to be the case. If accelerated elec- 
trons alone, unaccompanied by neutralizing ions, stream down the legs of the loop from the 
acceleration region to the footpoints, they will drive a co-spatial return current in the ambi- 
ent plasma to neutralize the high current associated with the downward- streaming electrons. 
We refer to this primarily downward-streaming distribution of energetic (suprathermal) elec- 
trons as an electron beam. The electric field associated with the return current will decelerate 
electrons in the beam, which can in turn modify the X-ray spectrum from the accelerated 
electrons. This is addressed in Sect. 5. 

Both the primary beam of accelerated electrons and the return current can become unsta- 
ble and drive the growth of waves in the ambient plasma. These waves can, in turn, interact 
with the electron beam and return current, altering the energy and angular distributions of 
the energetic electrons. These plasma instabilities are discussed in Sect. 6. 

The collisional energy loss rate is greater for lower energy electrons. Therefore, for 
suprathermal electrons streaming downward to the footpoints of a loop, the footpoint X- 
ray sources observed at lower energies should be at a higher altitude than footpoint sources 
observed at higher X-ray energies. The height dispersion of these sources provides infor- 
mation about the height distribution of the plasma density in the footpoints. The spatial 
resolution of the Reuven Ramaty High Energy Solar Spectroscopic Imager {RHESSI — see 
Lin et al. 2002) has made such a study possible. This is described in Sect. 7. RHESSI has 
observed X-ray sources that move downward from the loop top and then upward from the 
footpoints during some flares. This source evolution in time is also discussed in Sect. 7. 

If electrons of all energies are simultaneously injected, the footpoint X-ray emission 
from the slower, lower energy electrons should appear after that from faster, higher energy 
electrons. The length of this time delay provides an important test for the height of the accel- 
eration region. Longer time delays can result from magnetic trapping of the electrons. The 
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evolution of the thermal plasma in flares can also exhibit time delays associated with the bal- 
ance between heating and cooling processes. These various time delays and the information 
they provide are addressed in Sect. 8. 

An important diagnostic of electron acceleration and propagation in flares is the time 
evolution of the X-ray spectrum during flares. In most flares, the X-ray spectrum becomes 
harder (flatter, smaller spectral index) and then softer (steeper, larger spectral index) as the 
X-ray flux evolves from low to high intensity and then back to low intensity. There are 
notable exceptions to this pattern, however. Spectral evolution is addressed in Sect. 9. 

One of the most important results from the Yohkoh mission is the discovery in some 
flares of a hard (high energy) X-ray source above the top of the thermal (low energy) X-ray 
loops. This, together with the Yohkoh observations of cusps at the top of flare X-ray loops, 
provided strong evidence that energy release occurs in the corona above the hot X-ray loops 
(for some flares, at least). Although several models have been proposed, the origin of these 
“above-the-looptop” hard X-ray sources is not well understood. We need to determine how 
their properties and evolution compare to the more intense footpoint hard X-ray sources. 
These issues are addressed in Sect. 10. 

Radio observations provide another view of accelerated electrons and related flare phe- 
nomena. Although radio emission and its relationship to flare X-rays are primarily addressed 
in White et al. (201 1), some intriguing radio observations that bear upon electron accelera- 
tion in flares are presented in Sect. 1 1 . 

RHESSI observations of flare X-ray emission have led to substantial progress, but many 
questions remain unanswered. Part of the progress is that the questions are different from 
those that were asked less than a decade ago. The primary context for interpreting the X-ray 
emission from suprathermal electrons is still the thick-target model, but the ultimate goal is 
to understand how the electrons are accelerated. In Sect. 12 we summarize and discuss the 
implications of the X-ray observations for the thick-target model and electron acceleration 
mechanisms, and highlight some of the questions that remain to be answered. Implications 
of these questions for future flare studies are discussed. 

2 Thin- and Thick-Target X-ray Emission 

As was summarized in Sect. 1, the electron-ion bremsstrahlung X-rays from a beam of 
accelerated electrons will be most intense where the density of target ions is highest, as 
well as where the flux of accelerated electrons is high. The local emission (emissivity) 
at position r of photons of energy £ by electrons of energy E is given by the plasma 
ion density, w(r) ionscm“^, times the electron-beam flux density distribution, E(E,r) 
electrons cm“^ s“* keV^', times the differential electron-ion bremsstrahlung cross-section, 
2(fi, E) cm^ keV~*. For simplicity, we do not consider here the angular distribution of the 
beam electrons or of the emitted photons, topics addressed in Kontar et al. (2011). 

The emissivity of the radiation at energy e from all the electrons in the beam is obtained 
by integrating over all contributing electron energies, which is all electron energies above 
the photon energy. The photon flux emitted per unit energy is obtained by integrating over 
the emitting source volume (V) or, for an imaged source, along the line of sight through 
the source region. Finally, assuming isotropic emission, the observed spatially integrated 
flux density of photons of energy e at the X-ray detector, I (s) photons cm“^ s“* keV“', is 
simply the flux divided by the geometrical dilution factor 4jtR^, where R is the distance to 
the X-ray detector. Thus, 



(2.1) 
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We refer to 7(e) as the X-ray flux spectrum, or simply the X-ray spectrum. The spec- 
trum obtained directly from an X-ray detector is generally a spectrum of counts versus en- 
ergy loss in the detector, which must be converted to an X-ray flux spectrum by correcting 
for the detector response as a function of photon energy (see, for example. Smith et al. 
2002 ). 

Besides increasing the X-ray emission, a high plasma density also means increased 
Coulomb energy losses for the beam electrons. In the solar plasma, the bremsstrahlung 
losses are small compared to the collisional losses to the plasma electrons. For a fully ion- 
ized plasma and beam electron speeds much greater than the mean speed of the thermal 
electrons, the (nonrelativistic) energy loss rate is 


where K = 27re'* is the Coulomb logarithm for electron-electron collisions, e is the 

electron charge, «e(r) is the plasma electron number density, and v(7?) is the speed of the 
electron (see Brown 1971; Emslie 1978). The coefficient K is usually taken to be constant, 
although Aee depends weakly on the electron energy and plasma density or magnetic field 
strength, typically falling in the range 20-30 for X-ray-emitting electrons. Numerically, the 
energy loss rate in keV s“' or erg s“' , with E in keV is 


Here and in equations to follow, notation such as is used to show the scaling of 

a computed constant (here, X) with an independent variable and the numerical value taken 
for the independent variable (A^g = 23 in (2.3)). If the plasma is not fully ionized, K also 
depends on the ionization state (see Sect. 4.1). 

Noting that vr/t = dz, (2.2) can be simplified to dE/dNg = —KjE, where dNe(z) = 
rie{z)dz and Ne{z) (cm^^) is the plasma electron column density. (Here we treat this as a 
one-dimensional system and do not distinguish between the total electron velocity and the 
velocity component parallel to the magnetic field.) Hence, the evolution of an electron’s 
energy with column density is simply 


where Eq is the initial energy of the electron where it is injected into the target region. 
For example, a 1 keV electron loses all of its energy over a column density of 1/(2X) = 
1.7 X 10^^ cm“^ (for A^.^ = 23). A 25 keV electron loses 20% of its energy over a column 
density of 3.8 x 10*® cm“^. 

If energy losses are not significant within a spatially unresolved X-ray source region, 
the emission is called thin-target. If, on the other hand, the non-thermal electrons lose all 
their suprathermal energy within the spatially unresolved source during the observational 
integration time, the emission is called thick-target. We call a model that assumes these 
energy losses are from Coulomb collisions ((2.2) and (2.3)) a collisional thick-target model. 
Collisional thick-target models have been applied to flare X-ray/y-ray emission since the 
discovery of this emission in 1958 (Peterson and Winckler 1958, 1959). 

The maximum information that can be obtained about the accelerated electrons from an 
X-ray spectrum alone is contained in the mean electron flux distribution, the plasma-density- 
weighted, target-averaged electron flux density distribution (Brown et al. 2003; Kontar et al. 


dE/dt = -(K/E)nflr)\{E), 


(2.2) 



E^ = El~ 2KNe 


(2.4) 
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2011). The mean electron flux distribution is defined as 

F{E) = ^— f w(r) F(£', r) t/y electrons cm“^ keV^\ (2.5) 

nV Jy 

where n and V are the mean plasma density and volume of the emitting region. As can be 
seen from (2.1), the product hV F can, in principle, be deduced with only a knowledge of 
the bremsstrahlung cross-section, Q{e, E). Additional information is required to determine 
if the X-ray emission is thin-target, thick-target, or something in between. The flux distri- 
bution of the emitting electrons and the mean electron flux distribution are equivalent for a 
homogeneous, thin-target source region. 

Equation (2.1) gives the observed X-ray flux in terms of the accelerated electron flux 
density distribution throughout the source. However, we are interested in the electron dis- 
tribution injected into the source region, Fq(£'o, To), since that is the distribution produced 
by the unknown acceleration mechanism, including any modifications during propagation 
to the source region. To obtain this, we need to know how to relate E(E,r) at all locations 
within the source region to Tb(i?o, ro). Since we are interested in the X-rays from a spa- 
tially integrated, thick-target source region, the most direct approach is to first compute the 
bremsstrahlung photon yield from a single electron of energy Eq, v(s. Eg) photons keV^* 
per electron (Brown 1971). As long as the observational integration time is longer than the 
time required for the electrons to radiate all photons of energy e (i.e., longer than the time 
required for energy losses to reduce all electron energies to less than e), the thick-target 
X-ray spectrum is then given by 


^thick (^) 


1 

4nR^ 



.^o(Eo) v{e, Eo)dEo, 


( 2 . 6 ) 


where ^o(Eo) is the electron beam flux distribution (electrons s“' keV^*). ^o(Eo) is the 
integral of To(ro, Eo) over the area at the injection site through which the electrons stream 
into the thick-target region. 

The rate at which an electron of energy E radiates bremsstrahlung photons of energy e 
is n(r)v{E)Q{e, E). The photon yield is obtained by integrating this over time. Since the 
electrons are losing energy at the rate dE/dt, the time integration can be replaced by an 
integration over energy from the initial electron energy Eq to the lowest energy capable of 
radiating a photon of energy e: 


y(^, Eo) = 

Using (2.2) for dE/dt, (2.6) becomes 
1 1 


f 

JEq 


n{T)\(E)Q{e,E)dE 
dE/dt 


hhick )e) — 


ZK 


poo pEQ 

^ ^o{Eo) I EQ(e,E)dEdEo, 

J Eo=e J E=b 


(2.7) 


( 2 . 8 ) 


where we have used the relationship = Zn, where Z ~ 1.1 is the ion-species-number- 
density-weighted (or, equivalently, relative-ion-abundance-weighted) mean atomic number 
of the target plasma. 

By (2.8), the thick-target X-ray flux spectrum does not depend on the plasma density. 
However, the plasma must be dense enough for the emission to be thick-target, i.e., dense 
enough for all the electrons to be thermalized in the observation time interval. Integration of 
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(2.2) shows that this typically implies a plasma density > 10^*-10*^ cm“^ for an observa- 
tional integration time of 1 s (see Sects. 10.4 and 12.1 for more about this). This condition 
is well satisfied at loop footpoints. 

Observed non-thermal X-ray spectra from solar flares can usually be well fitted with a 
model photon spectrum that is either a single or a double power-law. For a single power-law 
electron flux distribution of the form ^(E) = AE~^, the photon spectrum is also well ap- 
proximated by the power-law form I (e) = IqS~'*' . The relationship between the electron and 
photon spectral indices S and y can most easily be obtained from (2.1) and (2.8) using the 
Kramers approximation or the nonrelativistic Bethe-Heitler (NRBH) bremsstrahlung cross 
section (see Koch and Motz 1959 for bremsstrahlung cross-sections). The NRBH cross- 
section is given by: 


Qnrbh{£, E) 


Z^Qo 

eE 


In 


1 + Vl -e/£^ 
1 - Vl - e/E ~ y 


cm^keV"*, 


(2.9) 


where Qo = 7.90 x 10“^^ cm^keV and ~ 1.4 is the ion-species-number-density- 
weighted mean square atomic number of the target plasma. The Kramers approximation 
to this cross-section is (2.9) without the logarithmic term. The bremsstrahlung cross-section 
is zero for e > £, since an electron cannot radiate a photon that is more energetic than 
the electron. Analytic expressions for the photon flux from both a uniform thin-target 
source and a thick-target source can be obtained with the Kramers and the NRBH cross- 
sections when the electron flux distribution has the single-power-law form (Brown 1971; 
Tandberg-Hanssen and Emslie 1988). The thin-target result also generalizes to the photon 
flux from a single-power-law mean electron flux distribution. 

For a uniform thin-target source and ^(E) = AE~^, 

I,Hin(s) = 3.93 X 10-^^ (2.10) 


giving Y,hi„ = 5 + 1 . The photon energy e is in keV and N is the ion column density. The 
power-law-index-dependent coefficient for the NRBH cross-section is 


P,n(S) 


B(S, 1/2) 
S 


( 2 . 11 ) 


where B(x, y) is the standard Beta function. (In the Kramers approximation, Ptn(S) = 1/5.) 
Typical values for the ion column density and A, the differential electron flux at 1 keV, are 
N = 10^*-10^° cm”^ and A = lO^'^-lO^^ electrons s“' keV“^ 

For a thick-target source region. 


W(£) = 1.17x 10-"" 





(2.12) 


giving Yihick = 5—1. The power-law-index-dependent coefficient for the NRBH cross- 
section is 


MS) 

In the Kramers approximation, Ptk(S) = 


_ B(S-2, 1/2) 

“ (5-I)(5-2)' 
l/[(5-I)(5-2)]. 


(2.13) 
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Fig. 1 The 

power-law-index-dependent 
terms fltn(S) (2.11) and Ptk(^) 

(2.13) in the analytic expressions 
for the thin-target (2.10) and the 
thick-target (2.12) photon flux 
from a power-law electron flux 
distribution. The solid curves are 
for the nonrelativistic 
Bethe-Heitler bremsstrahlung 
cross-section and the dotted 
curves are for the Kramers 
cross-section. Note that the fS 
axis is linear for the thin-target 
coefficient, and logarithmic for 
the thick-target coefficient 

2 4 6 8 10 12 

5 

The coefficients yS,„(i5) and Ptk{S) for both the NRBH and the Kramers cross-sections are 
plotted as a function of S in Fig. 1 . The Kramers and NRBH results are equal for thin-target 
emission when 5 3.4, and for thick-target emission when S ~ 5.4. For the plotted range 

of S, the Kramers approximation can differ from the NRBH result by over 90%. For S in 
the range 3-10, the Kramers result can differ from the NRBH result by as much as 76% and 
57% for thin- and thick-target emission, respectively. 

It is important to recognize that the above power-law relationships are only valid if the 
electron density distribution, F{r, E) electrons cm“^ s“* keV^' , or the electron dis- 
tribution, ^{E) electrons s“* keV“*, is assumed to have a power-law energy dependence. 
It is sometimes convenient to work with the electron density distribution, f(r,E) (elec- 
trons cm“^ keV“*), rather than the flux density distribution, especially when considering 
thin-target emission alone or comparing X-ray spectra with radio spectra. The flux density 
and density distributions are related through F{r, E) = /(r, E)v{E). If the electron density 
distribution rather than the flux or flux density distribution is assumed to have a power-law 
index S', so that f{r,E)(xE^^, the relationships between this power-law index and the 
photon spectral index become Yn,i„ = S' 0.5 and ythick = S' — 1.5. 

The simple power-law relationships are not valid if there is a break or a cutoff in the 
electron distribution at an energy less than ~2 orders of magnitude above the photon ener- 
gies of interest. Since all electrons with energies above a given photon energy e contribute 
to the bremsstrahlung at that photon energy, for the power-law relationships to be valid the 
break energy must be high enough that the deficit (or excess) of electrons above the break 
energy does not significantly affect the photon flux at energy e. The power-law relationship 
is typically not accurate until photon energies one to two orders of magnitude below the 
break energy, depending on the steepness of the power-law electron distribution (see Figs. 9 
and 10 of Holman 2003). Thus, for example, these relationships are not correct for the lower 
power-law index of a double power-law fit to a photon spectrum at photon energies within 
about an order of magnitude below the break energy in the double power-law electron dis- 
tribution. Equation (2.1) or (2.8) can be used to numerically compute the X-ray spectrum 
from an arbitrary flux distribution in electron energy. 

When electrons with kinetic energies approaching or exceeding 511 keV significantly 
contribute to the radiation, the relativistic Bethe-Heitler bremsstrahlung cross-section (Eq. 
3BN of Koch and Motz 1959) or a close approximation (Haug 1997) must be used. Haug 
(1997) has shown that the maximum error in the NRBH cross-section relative to the rela- 
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tivistic Bethe-Heitler cross section becomes greater than 10% at electron energies of 30 keV 
and above. Numerical computations using the relativistic Bethe-Heitler cross-section have 
been incorporated into the RHESSI spectral analysis software (OSPEX) for both thin- and 
thick-target emission from, in the most general case, a broken-power-law electron flux dis- 
tribution with both low- and high-energy cutoffs (the functions labeled “thin” and “thick” 
using the IDL programs brm_bremspec .pro and brm_bremthick.pro — see Hol- 
man 2003). Faster versions of these programs are now available in OSPEX and are currently 
labeled “thin2” and “thick2” and use the IDL programs brm2_thintarget .pro and 
brm2_thicktarget .pro. 

The analytic results based on the NRBH cross-section have been generalized to a broken- 
power-law electron flux distribution with cutoffs by Brown et al. (2008). They find a max- 
imum error of 35% relative to results obtained with the relativistic Bethe-Heitler cross- 
section for the range of parameters they consider. These results provide the fastest method 
for obtaining thin- and thick-target fits to X-ray spectra in the RHESSI spectral analysis soft- 
ware, where they are labeled “photon_thin” and “photon_thick” and use the IDL programs 
f photon thin. pro and f photon thick. pro. 

3 Low-Energy Cutoffs and the Energy in Non-thermal Electrons 

One of the most important aspects of the distribution of accelerated electrons is the low- 
energy cutoff. The acceleration of charged particles out of the thermal plasma typically 
involves a competition between the collisions that keeps the particles thermalized and the 
acceleration mechanism. The particles are accelerated out of the tail of the thermal distribu- 
tion, down to the lowest particle energy for which the acceleration mechanism can overcome 
the collisional force. Thus, the value of the low-energy cutoff can provide information about 
the force of the acceleration mechanism. More generally, as discussed below, the electron 
distribution must have a low-energy cutoff (1) so that the number and energy flux of elec- 
trons is finite and reasonable, and (2) because electrons with energies that are not well above 
the thermal energy of the plasma through which they propagate will be rapidly thermalized. 
Knowledge of the low-energy cutoff and its evolution during a flare is critical to determin- 
ing the energy flux and energy in non-thermal electrons and, ultimately, the efficiency of the 
acceleration process. 

3.1 Why Do We Need to Determine the Low-Energy Cutoff of Non-thermal Electron 
Distributions? 

An important feature of the basic thick-target model is that the photon spectrum 1(e) is di- 
rectly determined by the injected electron flux distribution ^o(Eo). As can be seen from 
(2.8), no additional parameters such as source density or volume need to be determined. 
Consequently, by integrating over all electron energies, we can also determine the total 
flux of non-thermal electrons, N„th electrons s“', the power in non-thermal electrons, Pnth 
erg s“* , and, integrating over time, the total number of, and energy in, non-thermal electrons. 

The total non-thermal electron number flux and power are computed as follows: 



■+00 


(3.1) 



-1 


(3.2) 


^ Springer 


G.D. Holman et al. 


Fig. 2 Typical full-Sun flare 
spectrum. Dashed: Nonthermal 
thick-target spectRim from an 
accelerated electron distribution 
with (5 = 4, and a low-energy 
cutoff of 20 keV. Dotted: 

Thermal spectrum, from a plasma 
with temperature T = 20 MK 
and emission measure 
EM = 10'*® cm“3. Solid: Total 
radiated spectrum. The multiple 
peaks in the thermal spectrum are 
from spectral lines, as observed 
by an instrument with ~ 1 keV 
spectral resolution 



Photon energy e [keV] 


The last expression in each equation is the result for a power-law electron flux distribution 
of the form ^q{Eq) = A ■ . The constant ke = 1.60 x 10“® is the conversion from keV 

to erg. Ec is a low-energy cutoff to the electron flux distribution. These expressions are 
valid and finite for S > 2 and Ec > 0. We call this form of low-energy cutoff a sharp low- 
energy cutoff. An electron distribution that continues below a transition energy Ec that has a 
positive slope, is flat, or in general has a spectral index Siou, < 1 also provides finite electron 
and energy fluxes, but these fluxes are somewhat higher than those associated with the sharp 
low-energy cutoff. 

For this single-power-law electron flux distribution with a sharp low-energy cutoff, the 
non-thermal power (erg s“*), and ultimately the non-thermal energy (erg), from the power- 
law electron flux distribution depends on only three parameters; 5, A, and Ec. Observations 
indicate that & is greater than 2 (Dennis 1985; Lin and Schwartz 1987; Winglee et al. 1991; 
Holman et al. 2003). Hence, were Ec = 0, the integral would yield an infinite value, a de- 
cidedly unphysical result! Therefore, the power-law electron distribution cannot extend all 
the way to zero energy with the same or steeper slope, and some form of low-energy cutoff 
in the accelerated electron spectrum must be present. As we will see, the determination of 
the energy at which this cutoff occurs is not a straightforward process, but it is the single 
most important parameter to determine (as the other two are generally more straightforward 
to determine — see Sect. 2 and Kontar et al. 2011). For example, with 8 = 4 (typical during 
the peak time of strong flares), a factor of 2 error in Ec yields a factor of 4 error in . For 
larger 5 (as found in small flares, or rise/decay phases of large flares), such an error quickly 
leads to an order of magnitude (or even greater) difference in the injected power P„,h and in 
the total energy in the non-thermal electrons accelerated during the flare I 

3.2 Why Is the Low-Energy Cutoff Difficult to Determine? 

The essence of the problem in many flare spectra is summarized in Fig. 2: the non-thermal 
power-law is well-observed above ~20 keV, but any revealing features that it might possess 
at lower energies, such as a low-energy cutoff, are masked by the thermal emission. 

Even if a spectrum does show a flattening at low energies that could be the result of a 
low-energy cutoff, other mechanisms that could produce the flattening must be ruled out (see 
Sect. 3.4). The low-energy cutoff has the characteristic feature, determined by the photon 
energy dependence of the bremsstrahlung cross-section (see (2.9)), that the X-ray spectrum 
eventually approaches a spectral index of y ~ 1 at low energies (cf. Holman 2003). It is 
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Fig. 3 Different shapes of low-energy cutoff in the injected electron distribution (left) lead to slightly dif- 
ferent photon spectra (right). The cutoff/tumover electron energy is Ec = 20 keV. The thin curve in the right 
panel demonstrates how the cutoff can be masked by emission from thermal plasma. See also Holman (2003) 
for a thorough discussion of bremsstrahlung spectra generated from electron power-laws with cutoff 


currently impossible, however, to observe a flare spectrum to low enough photon energies 
to see that it does indeed become this flat. Generally we can only hope to rule out the other 
mechanisms based on additional data and detailed spectral fits. 

3.3 What Is the Shape of the Low-Energy Cutoff, and How Does It Impact the Photon 
Spectrum and P^th ? 

Bremsstrahlung photon spectra are obtained from convolution integrals over the electron 
flux distribution ((2.1) and (2.8)). Hence, features in an electron distribution are smoothed 
out in the resulting photon spectrum (see, e.g.. Brown et al. 2006). 

As can be seen in Fig. 3, both a sharp cutoff at and a “turnover” (defined here to be 
a constant Fq{E) below Ec, a “plateau”) in the injected electron distribution lead to similar 
thick-target photon spectra. This subtle difference is difficult to discriminate observationally, 
and the problem is compounded by the dominance of the thermal component at low energies. 

A sharp cutoff would lead to plasma instabilities that should theoretically flatten the dis- 
tribution around and below the cutoff within microseconds (see Sect. 6). On the other hand, 
the electron flux distribution below the cutoff must be flatter than as demonstrated by 
(3.1), or the total electron number flux would be infinite. Having a constant value for the dis- 
tribution below Ec (turnover case) seems like a reasonable middle ground and approximates 
a quasilinearly relaxed electron distribution (Sect. 6; Krall and Trivelpiece 1973, Chap. 10). 
Coulomb collisional losses, on the other hand, yield an electron distribution that increases 
linearly at low energies (see Fig. 4), leading to a photon spectrum between the sharp cutoff 
case and the turnover case. 

Notice that the photon spectra actually flatten gradually to the spectral index of 1 at low 
energies from the spectral index of / = 5 -|- 1 at Fc and higher energies. Below Ec, it is not 
a power-law. Fitting a double power-law model photon spectrum, and using the break (i.e., 
kink) energy as the low-energy cutoff typically leads to a large error in Ec (e.g., Gan et al. 
2001; Saint-Hilaire and Benz 2005), and hence to an even larger error in P„th- 
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Fig. 4 The four plots show the Coulomb-collisional evolution with column density of an injected electron 
distribution {thick, solid line). For the simple power-law case (upper left), the low-energy end of the distri- 
bution becomes linear, and the peak of the distribution is found at = £* / sfE, where ^ is the injected 

distribution power-law spectral index (5 = 4 in the plots), and £* = ^2K ■ W* is the initial energy that elec- 
trons must possess in order not to be fully stopped by a column density Nn. (2.4). When a low-energy cutoff 
is present, the peak of the distribution is seen to first decrease in energy until £* exceeds the cutoff energy 
(from Saint-Hilaire 2005). 


In terms of the energetics, Saint-Hilaire and Benz (2005) have shown that the choice of 
an exact shape for the low-energy cutoff as a model is not dramatically important. For a fixed 
cutoff energy, from (3.2) it can be shown that the ratio of the power in the turnover model 
to the power in the sharp cutoff model without the flat component below the cutoff energy 
is 8/2. In obtaining spectral fits, however, the turnover model gives higher cutoff energies 
than the sharp cutoff model. Using simulations, Saint-Hilaire and Benz (2005) found that 
assuming either a sharp cutoff model or a turnover model led to differences in P„,h typi- 
cally less than ~20%. Hence, the sharp cutoff, being the simplest, is the model of choice 
for computing flare energetics. Nevertheless, knowing the shape of the low-energy cutoff 
would not only yield more accurate non-thermal energy estimates, but would be a source of 
information on the acceleration mechanism and/or propagation effects. 

Spectral inversion methods have recently been developed for deducing the mean electron 
flux distribution (2.5) from X-ray spectra (Johns and Lin 1992; Brown et al. 2003, 2006). A 
spectral “dip” has been found just above the presumed thermal component in some deduced 
mean electron flux distributions that may be associated with a low-energy cutoff (e.g., Piana 
et al. 2003). In the collisional thick-target model, the slope of the high-energy “wall” of this 
dip should be linear or flatter, with a linear slope indicating the absence of emitting electrons 
in the injected electron distribution at the energies displaying this slope. Kontar and Brown 
(2006a) have found evidence for slopes that are steeper than linear, but their spectra were not 
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corrected for photospheric albedo (see Sect. 3.5). Finding and understanding these dips is a 
crucial element for gaining an understanding of the low-energy properties of flare electron 
distributions (see Kontar et al. 2011). 

Emslie (2003) has pointed out that the non-thermal electron distribution could seam- 
lessly merge into the thermal distribution, removing the need for a low-energy cutoff. As 
was shown by Holman et al. (2003) for SOL2002-07-23T00:35 (X4.8), however, merger of 
the electron distribution into the typically derived ~ 10-30 MK thermal flare plasma gener- 
ally implies an exceptionally high energy in non-thermal electrons. Thus, for a more likely 
energy content, a higher low-energy cutoff or a hotter plasma would need to be present in 
the target region. Any emission from this additional “hot core,” because of its much lower 
emission measure, is likely to be masked by the usual ~ 10-30 MK thermal emission. This 
merger of the non-thermal electron distribution into the thermal tail in the target region does 
not remove the need for a low-energy cutoff in the electrons that escape the acceleration 
region, however. 

This section has dealt with the shape of the low-energy cutoff under the assumptions that 
the X-ray photon spectra are not altered by other mechanisms and that the bremsstrahlung 
emission is isotropic. The next section lists the important caveats to these assumptions, and 
their possible influence in the determination of the low-energy cutoff to the electron flux 
distribution. 

3.4 Important Caveats 

As previously discussed, apparently minor features in the bremsstrahlung photon spectrum 
can have substantial implications for the mean electron flux and, consequently, the injected 
electron distribution. This means that unknown or poorly-understood processes that alter 
the injected electron distribution (propagation effects, for example) or the photon spectrum 
(including instrumental effects) can lead to significant errors in the determination of the 
low-energy cutoff. Known processes that affect the determination of the low-energy cutoff 
are enumerated below. 

1 . Detector pulse pileup effects (Smith et al. 2002), if not properly corrected for, can in- 
troduce a flattening of the spectrum toward lower energies that simulates the flattening 
resulting from a low-energy cutoff. 

2. The contribution of Compton back-scattered photons (photospheric albedo) to the mea- 
sured X-ray spectrum can simulate the spectral flattening produced by a low-energy cut- 
off. Kasparova et al. (2005) have shown that the dip in a spectrum from SOL2002-08- 
20T08:25 (M3. 4) becomes statistically insignificant when the spectrum is corrected for 
photospheric albedo (also see Kontar et al. 2011). Kasparova et al. (2007) show that 
spectra in the 15-20 keV energy band tend to be flatter near disk center when albedo 
from isotropically emitted photons is not taken into account, further demonstrating the 
importance of correcting for photospheric albedo. 

3. The assumed differential cross-section and electron energy loss rate can influence the 
results (for a discussion of this, see Saint-Hilaire and Benz 2005). In some circumstances, 
a contribution from recombination radiation may significantly change the results (Brown 
et al. 2010; also see Kontar et al. 2011). 

4. Anisotropies in the electron beam directivity and the bremsstrahlung differential cross- 
section can significantly alter the X-ray spectrum (Massone et al. 2004). 

5. Non-uniform target ionization (the fact that the chromosphere’s ionization state varies 
with depth, see Sect. 4) can introduce a spectral break that may be confused with the 
break associated with a low-energy cutoff. 
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6. Energy losses associated with a return current produce a low-energy flattening of the X- 
ray spectrum (Sect. 5). This is a low-energy “cutoff” in the electron distribution injected 
into the thick target, but it is produced between the acceleration region and the emitting 
source region. 

7. A non-power-law distribution of injected electrons or significant evolution of the injected 
electron distribution during the observational integration time could affect the deduced 
value of the low-energy cutoff. 

For all the above reasons, the value of the low-energy cutoff in the injected electron flux 
distribution has not been determined with any degree of certainty except perhaps in a few 
special cases. Even less is known about the shape of the low-energy cutoff. The consensus 
in the solar physics community for now is to assume the simplest case, a sharp low-energy 
cutoff. Existing studies, presented in the next section, tend to support the adequacy of this 
assumption for the purposes of estimating the total power and energy in the accelerated 
electrons. 

3.5 Determinations of Ec and Electron Energy Content from Flare Data 

Before RHESSI, instruments did not cover well (if at all) the ~10^0 keV photon ener- 
gies where the transition from thermal emission to non-thermal emission usually occurs. 
Researchers typically assumed an arbitrary low-energy cutoff at a value at or below the in- 
strument’s observing range (one would talk of the “injected power in electrons above E^ 
keV” instead of the total non-thermal power Pnth)- An exception is Nitta et al. (1990). They 
argued that spectral flattening observed in two flares with the Solar Maximum Mission and 
Hinotori indicated a cutoff energy of >50 keV. Also, Gan et al. (2001) interpreted spec- 
tral breaks at ~80 keV in Compton Gamma-Ray Observatory (CGRO) flare spectra as the 
low-energy cutoff in estimating flare energetics, resulting in rather small values for the non- 
thermal energy in the analyzed flares. The relatively low resolution of the spectra from these 
instruments prevented the quantitative evaluation of any spectral flattening toward lower 
energies, however. 

The only high-resolution flare spectral data before the launch of RHESSI was the balloon 
data of Lin et al. (1981) for SOL1980-06-27T16:17 (M6.7) along with ~25 microflares ob- 
served during the same balloon flight. Benka and Holman (1994) applied a direct electric 
field electron acceleration model to the SOL 1980-06-27 flare data. They derived, along with 
other model-related parameters, the time evolution of the critical energy above which run- 
away acceleration occurs — the model equivalent to the low-energy cutoff. They found this 
critical energy to range from ~20-40 keV. 

It is now possible in most cases to obtain a meaningful upper limit on E^, thanks to 
RHESSI’s high-spectral-resolution coverage of the 10^0 keV energy range and beyond. 
Holman et al. (2003), Emslie et al. (2004), and Saint-Hilaire and Benz (2005), in determining 
the low-energy cutoff, obtained the “highest value for E^. that still fits the data.” In many solar 
flare spectra, because of the dominance of radiation from thermal plasma at low energies, 
a range of values for Ec fit the data equally well, up to a certain critical energy, above 
which the goodness-of-fit parameter becomes unacceptably large. The low-energy cutoff 
is taken to be equal to this critical value. This upper limit on the cutoff energy results in a 
lower limit for the non-thermal power and energy. The results obtained for the maximum 
value of Ec were typically in the 15^5 keV range, although late in the development of 
SOL2002-07-23T00:35 (X4.8) some values as high as ~80 keV were obtained for Ec- The 
minimum non-thermal energies thus determined were comparable to or somewhat larger 
than the calculated thermal energies. 


^ Springer 


Implications of X-ray Observations for Electron Acceleration 


1000.000 

100.000 

10.000 

1.000 

0.100 

0.010 

0.001 

1 



Photon Energy (keV) 


Photon Energy (keV) 


1000.000 

100.000 

10.000 

1.000 

0.100 

0.010 

0.001 

1 



Photon Energy (keV) 


Photon Energy (keV) 


Fig. 5 RHESSI spatially integrated spectra in four time intervals during SOL2002-04-15T03:55 (Ml. 2). 
(a) Spectrum at 23:06:20-23:06:40 UT (early rise phase), (b) Spectrum at 23:09:00-23:09:20 UT (just before 
impulsive phase), (c) Spectrum at 23:10:00-23:10:20 UT (soon after the impulsive rise), (d) Spectrum at 
23:11:00-23:11:20 UT (at the hard X-ray peak). The plus signs with error bars represent the spectral data. 
The lines represent model spectral fits: the dashed lines are non-thermal thick-target bremsstrahlung, the 
dotted lines are thermal bremsstrahlung, and the solid lines are the summation of the two (from Sui et al. 
2005a). Reproduced by permission of the AAS 


One of the best determinations of the low-energy cutoff so far was obtained by Sui et 
al. (2005a). They complemented the spatially-integrated spectral data for the SOL2002-04- 
15X03:55 (Ml. 2) limb flare with imaging and lightcurve information. Four spectra from 
this flare are shown in Fig. 5. The earliest spectrum, before the impulsive rise of the higher 
energy X-rays, was well fitted with an isothermal model. The last spectrum, from the time 
of the hard X-ray peak, clearly shows a thermal component below ~20 keV. Of particular 
interest is the second spectrum, showing both thermal and non-thermal fit components. As a 
consequence of the flattening of the isothermal component at low energies, the low-energy 
cutoff to the non-thermal component cannot extend to arbitrarily low energies without ex- 
ceeding the observed emission. This places a tight constraint on the value of the low-energy 
cutoff. The additional requirement that the time evolution of the derived temperature and 
emission measure of the thermal component be smooth and continuous throughout the flare 
constrains the value at other times. Applying the collisional thick-target model with a power- 
law distribution of injected electrons, they found the best cutoff value to be = 24 ± 2 keV 
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(roughly constant throughout the flare). The energy associated with these non-thermal elec- 
trons was found to be comparable to the peak energy in the X-ray-emitting thermal plasma, 
but an order of magnitude greater than the kinetic energy of the associated coronal mass ejec- 
tion (CME) (Sui et al. 2005b). This contrasts with results obtained for large flares, where the 
minimum energy in non-thermal electrons is typically found to be less than or on the order 
of the energy in the CME (e.g., Emslie et al. 2004). 

The importance of correcting for the distortion of spectra by albedo was revealed by a 
search for low-energy cutoffs in a sample of 177 flares with relatively flat spectra (y < 4) 
between 15 and 20 keV (Kontar et al. 2008a). Spectra can be significantly flattened by the 
presence of albedo photons in this energy range. The X-ray spectra, integrated over the 
duration of the impulsive phase of the flares, were inverted to obtain the corresponding 
mean electron flux distributions. Eighteen of the flares showed significant dips in the mean 
electron flux distribution in the 13-19 keV electron energy range that might be associated 
with a low-energy cutoff (see Sect. 3.3). However, when the X-ray spectra were corrected 
for albedo from isotropically emitted X-rays, all of the dips disappeared. Therefore, the 
authors concluded that none of these flare electron distributions had a low-energy cutoff 
above 12 keV, the lowest electron energy in their analysis. 

Low-energy cutoffs were identified in the spectra of a sample of early impulsive flares 
observed by RHESSI in 2002 (Sui et al. 2007). Early impulsive flares are flares in which 
the >25 keV hard X-ray flux increase is delayed by less than 30 s after the flux increase at 
lower energies. The pre-impulsive-phase heating of plasma to X-ray-emitting temperatures 
is minimal in these flares, allowing the nonthermal part of the spectrum to be observed 
to lower energies. In the sample of 33 flares, 9 showed spectral flattening at low energies 
in spectra obtained throughout the duration of each flare with a 4 s integration time. After 
correcting for the albedo from isotropically emitted X-rays, the flattening in 3 of the 9 flares, 
all near Sun center, disappeared. The flattening that persisted in the remaining 6 flares was 
consistent with that produced by a low-energy cutoff. The values derived for the low-energy 
cutoff ranged from 15 to 50 keV. The authors found the evolution of the spectral break and 
the corresponding low-energy cutoff in these flares to be correlated with the non-thermal 
hard X-ray flux. Further studies are needed to assess the significance of this correlation. 

Low-energy cutoffs with values exceeding 100 keV were identified in the spectra of 
the large flare SOL2005-01-19T08:22 (X1.3) (Warmuth et al. 2009). The hard X-ray light 
curve of this flare consisted of multiple peaks that have been interpreted as quasi-periodic 
oscillations driven by either magnetoacoustic oscillations in a nearby loop (Nakariakov et 
al. 2006) or by super-Alfvenic beams in the vicinity of the reconnection region (Ofman 
and Sui 2006). The high low-energy cutoffs were found in the last major peak of the series 
of hard X-ray peaks. Unlike the earlier peaks, this peak was also unusual in that it was not 
accompanied by the Neupert effect (see Sect. 8.3), consistent with the high values of the low- 
energy cutoff, and it exhibited soft-hard-harder rather than soft-hard-soft spectral evolution 
(see Sect. 9.1). A change in the character of the observed radio emission and movement 
of one of the two hard X-ray footpoints into a region of stronger photospheric magnetic 
field were also observed at the time of this peak. These changes suggest a strong connection 
between large-scale flare evolution and electron acceleration. 


4 Nonuniform Ionization in the Thick-Target Region 

In the interpretation of hard X-ray (HXR) spectra in terms of the thick-target model, one 
effect which has been largely ignored until recently is that of varying ionization along the 
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path of the thick-target beam. As first discussed by Brown (1973), the decrease of ionization 
with depth in the solar atmosphere reduces long-range collisional energy losses. This en- 
hances the HXR bremsstrahlung efficiency there, elevating the high energy end of the HXR 
spectrum by a factor of up to 2.8 above that for a fully ionized target. The net result is that 
a power-law electron spectrum of index S produces a photon spectrum of index y = S — I 
at low and high energies (see (2.12)), but with y < S — \ in between. The upward break, 
where the spectrum begins to flatten toward higher energies, occurs at fairly low energies, 
probably masked in measured spectra by the tail of the thermal component. The downward 
knee, where the spectrum steepens again to y = 3 — 1, occurs in the few deka-keV range, de- 
pending on the column depth of the transition zone. Thus, the measured X-ray spectrum may 
show a flattening similar to that expected for a low-energy cutoff in the electron distribution. 


4.1 Electron Energy Losses and X-ray Emission in a Nonuniformly Ionized Plasma 


The collisional energy-loss cross-section QdE) is dependent on the ionization of the back- 
ground medium. Flare-accelerated electron beams can propagate in the fully ionized corona 
as well as in the partially ionized transition region and chromosphere. Following Hayakawa 
and Kitao (1956) and Brown (1973), the cross-section QdE) can be written for a hydrogen 
plasma ionization fraction x 

'^1X6^ '2iJL 

QdE) = {xAee + (1 - x)Aen) = -^^A{x + X), (4.1) 

where e is the electronic charge, A^g the electron-electron logarithm for fully ionized media 
and Ach is an effective Coulomb logarithm for electron-hydrogen atom collisions. Numeri- 
cally Age ~ 20 and Agu ~ 7.1, so A = A^e — A^h — 12.9 and X = A^h! A Ci; 0.55. 

Then, in a hydrogen target of ionization level x(A) at column density A^(z), the energy 
loss equation for electron energy E is (cf. (2.2)) 

dE Ine^A K' 

— = ^(^+^(M)) = - (a + x(A)), (4.2) 

dN E E 

where K' = Ine^A = {AjA^dE ~ 0.65X. 

The energy loss of a given particle with initial energy Eq depends on the column density 
A!{z) = fg n{z')dz', so the electron energy at a given distance z from the injection site can 
be written E^ = E^ — 2K'M(N(z)) (cf. (2.4)), where 


N(z) 

M(N(z))= j {X + x(N'))dN' (4.3) 

0 

is the “effective” ionization-weighted collisional column density. 

The fractional atmospheric ionization x as a function of column density N (cm^^) 
changes from 1 to near 0 over a small spatial range in the solar atmosphere. Therefore, 
to lowest order, x(N) can be approximated by a step function x{N) = 1 for A < A*, and 
x{N) = 0 for A > A*. This gives M(N) = (A -|- 1)A for A < A* and M(N) = N^ + XN for 
A > A*. Electrons injected into the target with energies less than £* = y/2K'{X + 1)A* = 
a/2XA* experience energy losses and emit X-rays in the fully ionized plasma with x = I, 
as in the standard thick-target model. Electrons injected with energies higher than £* lose 
part of their energy and partially emit X-rays in the un-ionized (x = 0), or, more generally, 
partially ionized plasma. 
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Fig. 6 Photon spectrum 
residuals, normalized by the 
statistical eiTor for the spectral 
fit, for the time interval 
00:30:00-00:30:20 UT, 

2002-July-23, for {upper panel) 
an isothermal Maxwellian plus a 
power-law and {lower panel) an 
isothermal Maxwellian plus the 
nonuniform ionization spectrum 
with (5 = 4.24 and £* = 53 keV 
(from Kontar et al. 2003). 

Reproduced by permission of the 
AAS 

We can deduce the properties of the X-ray spectrum by substituting (4.2) into (2.7) (with 
dN = nvdt) and comparing Ithicki^) from (2.6) with Ithick(^) from (2.8). We see that for the 
nonuniformly ionized case the denominator in the inner integral now contains k + x{N) and 
K is replaced with K' . In the step-function model for j::(A^), photon energies greater than or 
equal to e* = ii* are emitted by electrons in the un-ionized plasma with £ > £*. Since k + 
x{N) has the constant value k, the thick-target power-law spectrum is obtained (for injected 
power-law spectrum), but the numerical coefficient contains K'k = Ine^Agn instead of K. 
At photon energies far enough below e* that the contribution from electrons with E > Etf 
is negligible, k + x{N) = 7. -|- 1 and the numerical coefficient contains {k+ \)K' = K. The 
usual thick-target spectral shape and numerical coefficient are recovered. The ratio of the 
amplitude of the high-energy power-law spectrum to the low-energy power-law spectrum is 
(k -b 1)/A. ~ 2.8. The photon energy e*(keV) ~ 2.3 x 10“^yN*(cm^2), between where the 
photon spectrum flattens below the high-energy power law and above the low-energy power 
law, determines the value of the column density where the plasma ionization fraction drops 
from 1 to 0. 

4.2 Application to Flare X-ray Spectra 

The step-function nonuniform ionization model was used by Kontar et al. (2002, 2003) to 
fit photon spectra from five flares. They assume a single power-law distribution of injected 
electrons with power-law index 8 and approximate the bremsstrahlung cross-section with 
the Kramers cross-section. First, they fit the spectra to the sum of a thermal Maxwellian 
at a single temperature T plus a single power law of index y ■ For SOL2002-07-23T00:35 
(X4.8) (Kontar et al. 2003) they limit themselves to deviations from a power law in the 
non-thermal component of the spectrum above ~40 keV. The top panel of Fig. 6 shows 
an example of such deviations, which represent significant deviations from the power-law 
fit. These deviations are much reduced by replacing the power law with the spectrum from 
the nonuniform ionization model, with the minimum rms residuals obtained for values of 
8 = 4.24 and £* = 53 keV (Fig. 6, bottom panel). The corresponding minimum (reduced) 

value obtained for the best fit to the X-ray spectrum (10-130 keV) dropped from 1.4 for 
the power-law fit to 0.8 for the nonuniform ionization fit. There are still significant residuals 
present in the range from 10 to 30 keV; these might be due to photospheric albedo or the 
assumption of a single-temperature thermal component. 

By assuming that the main spectral feature observed in a hard X-ray spectrum is due 
to the increased bremsstrahlung efficiency of the un-ionized chromosphere, allowance for 
nonuniform target ionization offers an elegant direct explanation for the shape of the ob- 
served hard X-ray spectrum and provides a measure of the location of the transition region. 



Energy (keV) 
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Table 1 Best fit nonuniformly ionized target model parameters for 
lent V* (energy range 20-100 keV), and the ratio of 

a single power-law the equiva- 

(from Kontar et al. 2002) 

Date 

Time, UT 

kT'(keV) 

S 

£* (keV) V* (cm2) 

■ 

^nonuni ‘ ^um 

20 Feb 2002 

11:06 

1.47 

5.29 

37.4 

2.7 X 1020 

0.032 

17 Mar 2002 

19:27 

1.27 

4.99 

24.4 

1.1 xlO^o 

0.047 

31 May 2002 

00:06 

2.02 

4.15 

56.2 

6.1 xlO^O 

0.041 

1 Jun 2002 

03:53 

1.45 

4.46 

21.0 

8.4 xlO'^ 

0.055 


Fig. 7 Variation of kT, S, £*, 
and Ntt throughout 
SOL2002-07-23T00:35 (X4.8) 
(Kontar et al. 2003). Reproduced 
by permission of the AAS. The 
variation of other parameters, 
such as emission measure, can be 
found in Holman et al. (2003) 
and Caspi and Lin (2010) 




Table 1 shows the best fit parameters derived for the four flare spectra analyzed by Kontar 
et al. (2002). The last column shows the ratio of the minimum value obtained from the 
nonuniform ionization fit to the minimum value obtained from a uniform ionization (sin- 
gle power-law) fit to the non-isothermal part of the spectrum. The nonuniform ionization 
model fits clearly provide substantially better fits than single power-law fits. 

Values of the fit parameters kT (keV), 8 and £* as a function of time for SOL2002- 
07-23T00:35 (X4.8), together with the corresponding value of lV*(cm"^) ~ 1.9 x 10'^£* 
(keV)^ were obtained by Kontar et al. (2003). The results (Fig. 7) demonstrate that the ther- 
mal plasma temperature rises quickly to a value ~3 keV and decreases fairly slowly there- 
after. The injected electron flux spectral index 8 follows a general “soft-hard-soft” trend and 
qualitatively agrees with the time history of the simple best-fit power-law index y (Holman 
et al. 2003). Z?* rises quickly during the first minute or so from ~40 keV to ~70 keV near 
the flare peak and thereafter declines rather slowly. The corresponding values of V* are 
~2-5 X 10^° cm-2. 

The essential results of these studies are that (1) for a single power-law electron injection 
spectrum, the expression for bremsstrahlung emission from a nonuniformly-ionized target 
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provides a significantly better fit to observed spectra than the expression for a uniform target; 
and (2) the value of f?* (and correspondingly A^*) varies with time. 

An upper limit on the degree of spectral flattening Ay that can result from nonuniform 
ionization was derived by Su et al. (2009). They applied this upper limit to spectra from a 
sample of 20 flares observed by RHESSI in the period 2002 through 2004. They found that 
15 of the 20 flare spectra required a downward spectral break at low energies and for each of 
these 15 spectra derived the difference Ay of the best-fit power-law spectral indices above 
and below the break. A Monte Carlo method was used to determine the 95% confidence 
interval for each of the derived values of Ay. Taking the value of Ay to be incompatible 
with nonuniform ionization if the 95% confidence interval fell above the derived upper limit, 
Su et al. (2009) found that six of the flare spectra could not be explained by nonuniform 
ionization alone. Thus, for these six flares some other cause such as a low-energy cutoff or 
return-current-associated energy losses (Sect. 5) must be at least partially responsible for the 
spectral flattening. 


5 Return Current Losses 

The thick-target model assumes that a beam of electrons is injected at the top of a loop 
and “precipitates” downwards in the solar atmosphere. Unless accompanied by an equal 
flux of positively charged particles, these electrons constitute a current and must cre- 
ate a significant self-induced electric field that in turn drives a co-spatial return cur- 
rent for compensation (Hoyng et al. 1976; Knight and Sturrock 1977; Emslie 1980; 
Diakonov and Somov 1988). The return current consists of ambient electrons, plus any pri- 
mary electrons that have scattered back into the upward direction. By this means we have a 
full electric circuit of precipitating and returning electrons that keeps the whole system neu- 
tral and the electron beam stable against being pinched off by the self-generated magnetic 
field required by Ampere’s law for an unneutralized beam current. However, the self-induced 
electric field results in a potential drop along the path of the electron beam that decelerates 
and, therefore, removes energy from the beam electrons. 

5.1 The Return Current Electric Field 

The initial formation of the beam/return-current system has been studied by van den Oord 
(1990) and references therein. We assume here that the system has time to reach a quasi- 
steady state. Van den Oord finds this time scale to be on the order of the thermal electron-ion 
collision time. This time scale is typically less than or much less than one second, depending 
on the temperature and density of the ambient plasma. In numerical simulations by Siversky 
and Zharkova (2009), times to reach a steady state after injection ranged from 0.07 s to 0.2 s, 
depending on the initial beam parameters. 

The self-induced electric field strength at a given location z along the beam and the flare 
loop, <o(z), is determined by the current density associated with the electron beam, j{z), 
and the local conductivity of the loop plasma, cr(z), through Ohm’s law: S’ (z) = j (z)/ct(z). 
Relating the current density to the density distribution function of the precipitating electrons, 
f (z, E,0), where E is the electron energy and 0 is the electron pitch angle, gives 



(5.1) 


0 0 
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Here /x is the cosine of the pitch angle and e and nij. are the electron charge and mass, 
respectively. The self-induced electric field strength S (z) depends on the local distribution 
of the beam electrons, which in turn depends on the electric field already experienced by 
the beam as well as any Coulomb energy losses and pitch-angle scattering that may have 
significantly altered the beam. It also depends on the local plasma temperature (and, to a 
lesser extent, density) through cr(z), which can, in turn, be altered by the interaction of the 
beam with the loop plasma (i.e., local heating and “chromospheric evaporation”). Therefore, 
determination of the self-induced electric field and its impact on the precipitating electrons 
generally requires self-consistent modeling of the coupled beam/plasma system. 

Such models have been computed by Zharkova and Gordovskyy (2005, 2006). They nu- 
merically integrate the time-dependent Fokker-Planck equation to obtain the self-induced 
electric field strength and electron distribution function along a model flare loop. The in- 
jected electron beam was assumed to have a single power-law energy distribution in the 
energy range from Eiow = 8 keV to E^pp = 384 keV and a normal (Gaussian) distribution 
in pitch-angle cosine ji with half-width dispersion A/x = 0.2 about /x = 1. 

The model computations show that the strength of the self-induced electric field is nearly 
constant at upper coronal levels and rapidly decreases with depth (column density) in the 
lower corona and transition region. The rapidity of the decrease depends on the beam 
flux spectral index. It is steeper for softer beams (S = 5-7) than for harder ones (5 = 3). 
The strength of the electric field is higher for a higher injected beam energy flux density 
(erg cm“^ s“ ' ), and the distance from the injection point over which the electric field strength 
is highest (and nearly constant) decreases with increasing beam flux density. 

5.2 Impact on Hard X-ray Spectra 

Deceleration of the precipitating beam by the electric field most significantly affects the 
lower energy electrons (<100 keV), since the fraction of the original particle energy lost 
to the electric field is greater for lower energy electrons. This leads to flattening of the 
electron distribution function towards the lower energies and, therefore, flattening of the 
photon spectrum. 

Photon spectra computed from kinetic solutions that include return current energy losses 
and collisional energy losses and scattering are shown in Fig. 8 (a) and (b). Low- and high- 
energy spectral indices and their dependence on the power-law index of the injected electron 
distribution and on the injected beam energy flux density are shown in Fig. 8 (c) and (d). 
The difference between the high-energy and low-energy spectral indices is seen to increase 
with both the beam energy flux density and the injected electron power-law index S. The 
low-energy index is found to be less than 2 for S as high as 5 when the energy flux density 
is as high as 10*^ erg cm“^ s“* . 

5.3 Observational Evidence for the Presence of the Return Current 

We have seen that return current energy losses can introduce curvature into a spectrum, 
possibly explaining the “break” often seen in observed flare X-ray spectra. A difficulty in 
directly testing this explanation is that the thick-target model provides the power (energy 
flux) in the electron beam (ergs“*), but not the energy flux density (ergcm'^ s“*). X-ray 
images provide information about the area of the target, but this is typically an upper limit 
on the area. Even if the source area does appear to be well determined, the electron beam 
can be filamented so that it does not All the entire area (the Ailing factor is less than 1). 
Also, if only an upper limit on the low-energy cutoff to the electron distribution is known. 
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Fig. 8 (a) Photon spectra computed from full kinetic solutions including return current losses and colli- 
sional losses and scattering. The top spectrum is for an injected single-power-law electron flux distribution 
between 8 keV and 384 keV with an index of 5 = 3, and the bottom spectrum is for S = 1. The injected 
electron energy flux density is 10^ ergcm”^ s“*. (b) Same as (a), but for an injected energy flux density of 
10*^ ergcm“^ s“* . The tangent lines at 20 and 100 keV demonstrate the determination of the low-energy 
and high-energy power-law spectral indices yjow and Ykigh ■ (c) The photon spectral indices yiow {dashed 
lines) and yhigh {solid lines) vs. S for an injected energy flux density of 10® {squares), lO*** {circles), and 
10*^ ergcm“^ s“' (crosses), (d) yhigh~Ylow the log of the injected electron energy flux density for S 
equal to 3 {bottom curve, squares), 5 {middle curve, circles), and 7 {top curve, triangles) (from Zharkova and 
Gordovskyy 2006). Reproduced by permission of the AAS 


as described in Sect. 3.5, the energy flux density may be higher. Therefore, the observations 
typically only give a lower limit on the beam energy flux density. 

The non-thermal hard X-ray flux is proportional to the electron beam flux density, but 
the return-current energy losses are also proportional to the beam flux density. As a conse- 
quence, Emslie (1980) concluded that the flux density of the non-thermal X-ray emission 
from a flare cannot exceed on the order of 10“^^ cm“^ s“^ above 20 keV. Alexander and 
Daou (2007) have deduced the photon flux density from non-thermal electrons in a sample 
of 10 flares ranging from GOES class Ml. 8 to X17. They find that the non-thermal photon 
flux density does not monotonically increase with the thermal energy flux, hut levels off (sat- 
urates) as the thermal energy flux becomes high. They argue that this saturation most likely 
results from the growing importance of return-current energy losses as the electron beam 
flux increases to high values in the larger flares. They find that the highest non-thermal 
photon flux densities agree with an upper limit computed by Emslie. 

A correlation between the X-ray flux and spectral break energy was found by Sui et al. 
(2007) in their study of X-ray spectra in early impulsive flares (see Sect. 3.5). They point 
out that the increasing impact of return current energy losses on higher energy electrons as 
the electron beam energy flux density increases could be an explanation for this correlation. 
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Battaglia and Benz (2008) studied two flares with non-thermal coronal hard X-ray 
sources for which the difference between the measured photon spectral index at the foot- 
points and the spectral index of the coronal source was greater than two, the value expected 
for coronal thin-target emission and footpoint thick-target emission from a single power-law 
electron distribution (see Sect. 10.2). They argue that return-current losses between the coro- 
nal and footpoint source regions are most likely responsible for the large difference between 
the spectral indices. 

The return current can also affect the spectral line emission from flares. Evidence for 
the presence of the return current at the chromospheric level from observations of the linear 
polarization of the hydrogen Ha and H/i lines has been presented by Henoux and Karlicky 
(2003). Dzifcakova and Karlicky (2008) have shown that the presence of a return current 
in the corona may have a distinguishable impact on the relative intensities of spectral lines 
emitted from the corona. 


6 Beam-Plasma and Current Instabilities 

Interaction of the accelerated electrons with plasma turbulence as they stream toward the 
thick-target emission region can modify the electron distribution. In this section we briefly 
discuss a likely source of such turbulence: that generated by the electron beam itself. If the 
beam is or becomes unstable to driving the growth of plasma waves, these waves can interact 
with the beam and modify its energy and/or pitch angle distribution until the instability 
is removed or the wave growth is stabilized. The return current associated with the beam 
(Sect. 5) can also become unstable, resulting in greater energy loss from the beam. Beam- 
plasma and return-current instabilities in solar flares have been reviewed by Melrose (1990) 
and Benz (2002). 

A sharp, low-energy cutoff to an electron beam or, more generally, a positive slope in 
the beam electron energy distribution is well-known to generate the growth of electrostatic 
plasma waves (Langmuir waves). The characteristic time scale for the growth of these waves 
is on the order of [(Mfc/w)o)pe]^*, where Wpe is the electron plasma frequency and Nb/n the 
ratio of the density of unstable electrons in the beam to the plasma density. This is on the 
order of microseconds for a typical coronal loop plasma density and ~ 10“^. This 
plasma instability is often referred to as the bump-on-tail instability. The result is that on 
a somewhat longer but comparable time scale electrons from the unstable region of the 
electron distribution lose energy to the waves until the sharp cutoff is flattened so that the 
distribution no longer drives the rapid growth of the waves. Therefore, the electron energy 
distribution below the low-energy cutoff is likely to rapidly become flat or nearly flat (suffi- 
ciently flat to stabilize the instability) after the electrons escape the acceleration region (see 
Chaps. 9 and 10 of Krall and Trivelpiece 1973). 

A recent simulation of the bump-on-tail instability for solar flare conditions, including 
Coulomb collisions and wave damping as the electrons propagate into an increasingly dense 
plasma, has been carried out by Hannah et al. (2009). The authors compute the mean electron 
flux distribution in their model flare atmosphere after injecting a power-law distribution 
with a sharp low-energy cutoff. They find a mean electron flux distribution with no dip (see 
Sect. 3.3) and a slightly negative slope below the original cutoff energy with a spectral index 
8 between 0 and 1 . 

A beam for which the mean electron velocity parallel to the magnetic held substantially 
exceeds the mean perpendicular velocity can drive the growth of waves that resonantly in- 
teract with the beam. When the electron gyrofrequency exceeds the plasma frequency, these 
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waves are electrostatic and primarily scatter the electrons in pitch angle. Generally known 
as the anomalous Doppler resonance instability, this instability tends to isotropize the beam 
electrons. Holman et al. (1982) showed that under solar flare conditions this instability can 
grow and rapidly isotropize the beam electrons in less than a millisecond. They found that 
electrons at both the low- and high-energy ends of the distribution may remain unscattered, 
however, because of wave damping. This could result in up to two breaks in the emitted X- 
ray spectrum. On the other hand, Vlahos and Rowland (1984) have argued that non-thermal 
tails will form in the ambient plasma and stabilize the anomalous Doppler resonance insta- 
bility by suppressing the growth of the plasma waves. 

Electrons streaming into a converging magnetic field can develop a loss-cone distribu- 
tion, with a deficit of electrons at small pitch angles. Both classical Coulomb collisions and 
loss-cone instabilities can relax this distribution by scattering electrons into the loss cone or 
extracting energy from the component of the electron velocities perpendicular to the mag- 
netic field (e.g., Aschwanden 1990). One loss-cone instability, the electron-cyclotron or gy- 
rosynchrotron maser, produces coherent radiation observable at radio frequencies (Holman 
et al. 1980; Melrose and Dulk 1982). 

The return current associated with the streaming electrons becomes unstable to the ion- 
acoustic instability when its drift speed exceeds a value on the order of the ion sound speed. 
The excited ion sound waves enhance the plasma resistivity, increasing the electric field 
strength associated with the return current, the heating of the plasma by the current, and the 
energy loss from the electron beam. 

It has been argued that rapid plasma heating and particle acceleration in the corona should 
result in the expansion of hot plasma down the legs of flare loops at the ion sound speed, 
confined behind a collisionless ion-acoustic conduction front (Brown et al. 1979). Electrons 
with speeds greater than about three times the electron thermal speed would be able to stream 
ahead of the conduction front. This scenario has not been observationally verified, but the 
observational signature may be confused by the chromospheric evaporation produced by the 
high-energy particles streaming ahead of the conduction front. 

Rowland and Vlahos (1985) argued that if the electron beam is unstable to beam plasma 
interactions, the return current will be carried by high-velocity electrons. This reduces the 
impact of collisions on the beam/return-current system and helps stabilize the system. In a 
recent simulation, Karlicky et al. (2008) have found that for current drift velocities exceeding 
the electron thermal speed, the return current is carried by both the primary (drifting thermal) 
current and an extended tail of high-velocity electrons. 

The evolution of the electron-beam/retum-current system when the return-current drift 
speed exceeds the electron thermal speed has also been simulated by Lee et al. (2008). 
They find that double layers form in the return current, regions of enhanced electric field 
that further increase the energy losses of the electron beam. This, in turn, increases the 
highest electron energy to which these losses significantly flatten the electron distribution 
and corresponding hard X-ray spectrum. 

The beam/return-current system has been simulated by Karlicky (2009), with a focus on 
the role of the Weibel instability. The Weibel instability tends to isotropize the electron distri- 
bution. Karlicky and Kasparova (2009) have computed the thin-target X-ray emission from 
the evolved electron distributions for a model with a weak magnetic field and another model 
with a strong magnetic field (ratio of the electron gyrofrequency to the plasma frequency 
~0 and ~1, respectively). They demonstrate that in both cases the electron distribution is 
more isotropic and the directivity of the X-ray emission is lower than when the instability of 
the system is not taken into account, with the greatest isotropization occurring in the weak 
field limit. 


^ Springer 


Implications of X-ray Observations for Electron Acceleration 


Although we expect plasma instabilities to affect the evolution of the electron beam, ob- 
servationally identifying them is difficult. The bump-on-tail instability and return current 
losses both lead to a flat low-energy cutoff. So far we have not established the ability to 
observationally distinguish a flat low-energy cutoff from a sharp low-energy cutoff. The 
bump-on-tail instability may be distinguishable from return-current losses by its short time 
scale and, therefore, the short distance from the acceleration region over which it effectively 
removes the unstable positive slope from the electron energy distribution. The instabilities 
that isotropize the electron pitch-angle distribution may be responsible for evidence from 
albedo measurements that flare electron distributions are isotropic or nearly isotropic (Kon- 
tar and Brown 2006b). 


7 Height Dependence and Size of X-ray Sources with Energy and Time 

7.1 Footpoint Sources 

Hard X-ray footpoint sources result from collisional bremsstrahlung radiated by precipitat- 
ing electrons, which produce most of the emission in the chromosphere according to the 
collisional thick-target model. Depending on the density structure in the legs of the coronal 
magnetic loop, mildly energetic electrons lose their energy in the lower corona or transi- 
tion region, while the more energetic electrons penetrate deeper into the chromosphere (see 
(2.2)). 

The altitude of these hard X-ray footpoint sources could never be measured accurately 
before RHESSI, because of a lack of spatial and spectral resolution. With RHESSI, we can 
measure the centroid of the footpoint location with an accuracy of order an arcsecond for 
every photon energy in steps as small as 1 keV. For a flare near the limb (Fig. 9), the centroid 
location translates directly into an altitude. 

Aschwanden et al. (2002) studied such a flare, SOL2002-02-20T1 1:07 (C7.5). The 
heights of the footpoint sources were fitted with a power-law function of the photon energy. 
This yielded altitudes h ~ 1000-5000 km in the energy range e = 10-60 keV, progressively 
lower with higher energy, as expected from the thick-target model (Fig. 9, right frame). 

Since the stopping depth of the precipitating electrons is a function of column density, 
the integrated density along their path in the chromosphere (2.4), the measured height de- 
pendence of the hard X-ray centroids can be inverted to yield a density model of the chromo- 
sphere (Brown et al. 2002). Assuming the decrease in density with height had a power-law 
dependence and the plasma is fully ionized, the inversion of the RHESSI data in the example 
shown in Fig. 9 yielded a chromospheric density model that has a significantly higher elec- 
tron density in the h = 2000-5000 km range than the standard chromospheric models based 
on UV spectroscopy and hydrostatic equilibrium (VAL and FAL models). The RHESSI- 
based chromospheric density model was therefore found to be more consistent with the 
“spicular extended chromosphere,” similar to the results from sub-mm radio observations 
during solar eclipses carried out at Caltech (Ewell et al. 1993). 

Forward fitting RHESSI X-ray visibilities to an assumed circular Gaussian source shape, 
Kontar et al. (2008b) found for a limb flare the full width at half maximum (FWHM) size 
and centroid positions of hard X-ray sources as a function of photon energy with a claimed 
resolution of ~0".2. They show that the height variation of the chromospheric density and 
of the magnetic flux density can be found with a vertical resolution of ~ 150 km by mapping 
the 18-250 keV X-ray emission of energetic electrons propagating in the loop at chromo- 
spheric heights of 400-1500 km. Assuming collisional losses in neutral hydrogen with an 
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Fig. 9 The centroids of footpoint hard X-ray emission are marked for different photon energies between 
10 keV and 60 keV for SOL2002-02-20T11:07 (C7.5), which occurred near the solar west limb and was 
imaged with RHESSI (left panel). The altitude /j(e) as a function of energy 8 shows a systematic height 
decrease with increasing energy (right panel) (from Aschwanden et al. 2002) 


exponential decrease in density with height, their observations of SOL2004-0 1-06X06:29 
(M5.8) suggest that the density of the neutral gas is in good agreement with hydrostatic 
models with a scale height of around 140 ± 30 km. FWHM sizes of the X-ray sources de- 
crease with energy, suggesting the expansion (fanning out) of magnetic flux tubes in the 
chromosphere with height. The magnetic scale height B(z){dB /dz)~^ is found to be on the 
order of 300 km and a strong horizontal magnetic field is associated with noticeable flux 
tube expansion at a height of ~900 km. A subsequent analysis with an assumed elliptical 
Gaussian source shape (Kontar et al. 2010) confirms these results and shows that the verti- 
cal extent of the X-ray source decrease with increasing X-ray energy. The authors find the 
vertical source sizes to be larger than expected from the thick-target model and suggest that 
a multi-threaded density structure in the chromosphere is required. The thick-target model 
to which the results were compared, however, did not account for partial occultation of the 
X-ray sources by the solar limb. 

The flare SOL2002-02-20T1 1 :07 (C7.5) has been reanalyzed by Prato et al. (2009) using 
both photon maps over a range of photon energies and mean electron flux maps deduced 
from RHESSI visibilities over a range of electron energies. Using source centroids computed 
from the maps and assuming an exponential decrease in density with height, they found the 
density scale height to be an order of magnitude larger than the expected chromospheric 
scale height on the quiet Sun, but consistent with the scale height in a non-static, flaring 
atmosphere. This is also consistent with the enhanced plasma densities found at ~1000- 
5000 km altitudes by Aschwanden et al. (2002). 

If the results for the 400-1500 km height range (Kontar et al. 2008b) and for the ~1000- 
5000 km height range (Aschwanden et al. 2002; Prato et al. 2009) are typical of flare 
loops, they imply that the upper chromosphere and transition region respond with a non- 
hydrostatic, expanded atmosphere while the low chromosphere does not respond to the flare 
energy release. These results could, of course, depend on the magnitude of the flare. More 
studies of this kind are clearly desirable, especially in coordination with observations of 
spectral lines from the chromosphere and transition region. 
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Fig. 10 CLEAN images at 
04:58:22-04:58:26 UT during the 
impulsive phase of 
SOL2003-11-13T05:01 (MI.6). 
The background shows the image 
at 9-12 keV. The contour levels 
are 75% and 90% of the peak 
flux at 9-12 keV {looptop), 70% 
and 90% at 12-18 keV (legs), 
and 50%, 60%, & 80% at 
28^3 keV ifootpoints) (from Liu 
et al. 2006). Reproduced by 
permission of the AAS 
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7.2 Loop Sources and Their Evolution 

As discussed above, footpoint sources are produced by bremsstrahlung emission in the thick- 
target chromosphere. The compactness of such sources results from the rapid increase of the 
density from the tenuous corona to the much denser chromosphere. This also gives rise 
to the compact height distribution of emission centroids at different energies as shown in 
Fig. 9. However, if the density distribution has a somewhat gradual variation, one would 
expect a more diffuse height distribution. Specifically, at some intermediate energies, we 
expect that HXR emission would appear in the legs of the loop, rather than the commonly 
observed looptop sources at low energies and footpoint sources at high energies. This has 
been observed by RHESSI in SOL2003- 1 1 - 1 3T05 :0 1 (M 1 .6) (Liu et al. 2006) (Fig. 1 0) and 
in SOL2002-1 1-28T04:37 (Cl.O) (Sui et al. 2006b). 

To reveal more details of the energy-dependent structure of SOL2003-11-13T05:01 
(M1.6), Figs, lla-c show the X-ray emission profile along the flare loop at different en- 
ergies for three time intervals in sequence. The high energy emission is dominated by the 
footpoints, but there is a decrease of the separation of the footpoints with decreasing energy 
and with time. At later times the profile becomes a single source, peaking at the looptop. 
The general trend suggests an increase of the plasma density in the loop with time (Liu et 
al. 2006), which can be produced by chromospheric evaporation and can give rise to pro- 
gressively shorter stopping distances for electrons at a given energy. Such a density increase 
also smooths out to some extent the sharp density jump at the transition region. This results 
in the non-thermal bremsstrahlung HXRs at intermediate energies appearing in the legs of 
the loop, at higher altitudes than the footpoints, as shown in Fig. 10. 

From the emission profiles in the non-thermal regimes of the photon spectra, Liu et 
al. (2006) derived the density distribution along the loop, using the empirical formula for 
non-thermal bremsstrahlung emission profiles given by Leach and Petrosian (1983, Eq. 1 1). 
Leach and Petrosian found that this formula closely approximates their numerical results 
for a steady-state, power-law injected electron distribution with a uniform pitch-angle dis- 
tribution, no return-current losses, and a loop with no magnetic field convergence. Since this 
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Fig. 11 (a) Brightness profiles in different energy bands measured along a semi-circular path fit to the flaring 
loop for the time interval 04:58:00-04:58:24 UT of SOL2003-11-13T05:01. The vertical axis indicates the 
average photon energy (logarithmic scale) of the energy band for the profile. Representative energy bands 
(in units of keV) are labeled above the corresponding profiles. The filled circles mark the local maxima, and 
the vertical dotted lines are the average positions of the centroids of the looptop and footpoint sources, (b, c) 
Same as (a), but for 04:58:24-04:58:48 and 04:58:48-04:59:12 UT, respectively. The error bars show the 
uncertainty of the corresponding profile (from Liu et al. 2006). Reproduced by permission of the AAS 


^ Springer 


Implications of X-ray Observations for Electron Acceleration 


Fig. 12 Averaged density 
profiles along one leg of the loop 
infen'ed from the HXR brightness 
profiles during the three time 
intervals in Fig. 1 1 . The distance 
is measured along the leg 
extending from the centroid of 
the thermal looptop source at 
about 15 arcsec in Fig. 1 1 to the 
end of the fitted semi-circle at 
about 37 arcsec (from Liu et al. 
2006). Reproduced by 
permission of the AAS 
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formula is a function of the column density, one does not need to assume any model form of 
the density distribution (cf. Aschwanden et al. 2002). Figure 12 shows the density profiles 
derived from the emission profiles in the three time intervals shown in Fig. 1 1 . Between the 
first and second intervals, the density increases dramatically in the lower part of the loop, 
while the density near the looptop remains essentially unchanged. The density enhance- 
ment then shifts to the looptop from the second to the third interval. This indicates a mass 
flow from the chromosphere to the looptop, most likely caused by chromospheric evapora- 
tion. For papers studying chromospheric evaporation using coordinated RHESSI FIXR and 
EUV Doppler-shift observations, see Milligan et al. (2006a, 2006b) and Brosius and Holman 
(2007). 

The flare SOL2002-1 1-28T04:37 (Cl.O) was an early impulsive flare, meaning that there 
was minimal pre-heating of plasma to X-ray-emitting temperatures prior to the appearance 
of impulsive hard X-ray emission (see Sect. 3.5). RHESSI observations of this flare showed 
coronal X-ray sources that first moved downward and then upward along the legs of the 
flare loop (Sui et al. 2006b). The bottom panel of Fig. 13 shows the motion of the sources 
observed in the 3-6 keV band. RHESSI and GOES light curves are shown in the top panel. 
The sources originated at the top of the flare loop and then moved downward along both legs 
of the loop until the time of peak emission at energies above 12 keV. Afterward the source 
in the northern leg of the loop was no longer observable, but the source in the southern 
leg moved back to the top of the loop. Its centroid location at the looptop was slightly but 
significantly lower than the centroid position at the beginning of the flare. Higher-energy 
sources showed a similar evolution, but they had lower centroid positions than their lower 
energy counterparts, again in agreement with the predictions of the thick-target model. 

The early downward source motion along the legs of the loop is a previously unobserved 
phenomenon. At this time we do not know if the occurrence is rare, or if it is simply rarely 
observed because of masking by the radiation from the thermal plasma. Sui et al. (2006b) 
argue that the motion results from the hardening of the X-ray spectrum, and possibly an 
increase in the low-energy cutoff, as the flare hard X-ray emission rises to its peak intensity. 
A flatter spectrum results in a higher mean energy of the electrons contributing to the radia- 
tion at a given X-ray energy. In a loop with a plasma density that increases significantly from 
the top to the footpoints, these higher energy electrons will propagate to a lower altitude in 
the loop as the spectrum hardens. The softening of the spectrum after peak emission would 
also contribute to the upward motion of the source after the peak. However, at that time chro- 
mospheric evaporation would likely be increasing the density in the loop, as discussed above 
for SOL2003-11-13T05:01 (Ml. 6), and thermal emission would be more important. All of 
these can contribute to an increase in the height of the centroid of the X-ray source. The 
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Fig. 13 RHESSI {solid lines) and GOES 1-8 A (dotted line) light curves are shown in the top panel. The 
RHESSI energy bands (from top to bottom) are 3-6, 6-12, 12-25, and 50-100 keV, with scaling factors of 
5, 1, 4, 3, and 0.5, respectively. The RHESSI and GOES integration times are 4 and 3 s, respectively. The 
bottom panel shows the distance between the 3-6 keV moving source centroids and their corresponding 
footpoint centroids located in the 25-50 keV image of the flare at the time of peak emission. The distances 
are plane-of-sky values with no correction for motions away from or toward the observer (from Sui et al. 
2006b). Reproduced by permission of the AAS 


downward motion may only occur in initially cool flare loops, i.e., early impulsive flares, 
because these loops are most likely to contain the density gradients that are required. 

In an attempt to differentiate between thermal and non-thermal X-ray emission, Xu et al. 
(2008) modeled the size dependence with photon energy of coronal X-ray sources observed 
by RHESSI in ten M-class limb flares. They determined the one-sigma Gaussian width of the 
sources along the length of the flare loops by obtaining forward fits to the source visibilities. 
The integration times ranged from one to ten minutes and the source sizes were determined 
in up to eight energy bins ranging in energy from as low as 7 keV to as high as 30 keV. They 
found the source sizes to increase slowly with photon energy, on average as The results 
were compared with several models for the variation of the source size with energy. The 
source size was expected to vary as for a thermal model with a constant loop density 
and a temperature that decreased with a Gaussian profile along the legs of the loop from 
a maximum temperature at the top of the loop. For the injection of a power-law electron 
flux distribution into a high-density loop so that the loop is a collisional thick target, the 
source size was expected to increase as e^. Neither of these models are consistent with the 
observed dependence. A hybrid thermal/non-thermal model and a non-thermal model 
with an extended acceleration region at the top of the loop were found to be consistent 
with the deduced scaling, however. The extended acceleration region was deduced to have 
a half-length in the range 10"-18" and density in the range (1-5) x 10" cm“^. We note 
that the extended acceleration region model implies a column density in the range 0.73- 
6.5 X 10^° cm“^ along the half length and, from (2.4), all electrons with energies less than 
somewhere in the range of 23-68 keV that traverse this half length will lose all of their 
energy to collisions. The acceleration process would therefore need to be efficient enough to 
overcome these losses. On the other hand, the 7-30 keV energy range is the range in which 
fits to spatially integrated X-ray spectra typically show a combination of both thermal and 
non-thermal bremsstrahlung emission. 
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Studies of flare hard X-ray source positions and sizes as a function of photon energy and 
time hold great promise for determining the height structure of flare plasma and its evolution, 
as well as providing information about the magnetic structure of the flare loop. Such studies 
are currently in their early stages, in that they usually assume an over simplified power-law 
or exponential height distribution for the plasma and do not take into account the variation 
of the plasma ionization state with height. They also assume the simple, one-dimensional 
collisional thick-target model, without consideration of the pitch-angle distribution of the 
beam electrons or the possibility of additional energy losses to the beam (such as return- 
current losses). Given the potential for obtaining a better understanding of flare evolution, 
we look forward to the application of more sophisticated models to the flare hard X-ray data. 


8 Hard X-ray Timing 

The analysis of energy-dependent time delays allows us to test theoretical models of physi- 
cal time scales and their scaling laws with energy. In the wavelength domain of hard X-rays, 
there are at least three physical processes known in the observation of solar flares that lead 
to measurable time delays as a function of energy (for a review, see Aschwanden 2004): 
(1) time-of-flight dispersion of free-streaming electrons, (2) magnetic trapping with the col- 
lisional precipitation of electrons, and (3) cooling of the thermal plasma. 

8.1 Time-of-Flight Delays 

The first type, the time-of-flight (TOF) delay, has a scaling of Af (e) oc and is caused 
by velocity differences of electrons that propagate from the coronal acceleration site to the 
chromospheric energy-loss region. The time differences are of order At ~ 10-100 ms for 
non-thermal electrons at energies E ~ 20-100 keV (e.g., Aschwanden et al. 1995, 1996). 
The measurement of such tiny time delays requires high photon statistics and high time 
resolution. Such data were provided by CGf?0/BATSE, which had 8 detectors, each with an 
effective collecting area of ~2000 cm^ and oriented at different angles to the Sun so that 
detector saturation at high count rates was not a problem. (For comparison, the total effective 
collecting area of RHESSI’s detectors is less than 100 cm^.) 

These studies of TOF delays have provided important evidence that electrons are accel- 
erated in the corona, above the top of the hot flare loops observed in soft X-rays. The fine 
structure in the light curves of most, but not all, of the studied flare bursts showed energy- 
dependent time delays consistent with the free streaming of electrons to the footpoints of 
the flare loops from an origin somewhat more distant than the half-length of the loops (As- 
chwanden et al. 1995, 1996; Aschwanden and Schwartz 1995). 

8.2 Trapping Delays 

The second type, the trapping delay, is caused by magnetic mirroring of coronal electrons 
which precipitate toward the chromosphere after a collisional time scale At(e) oc This 
is observable for time differences of At ~ 1-10 s for non-thermal electrons at E ^ 20- 
100 keV (e.g., Vilmer et al. 1982; Aschwanden et al. 1997). For trapping delays the higher 
energy X-rays lag the lower energy X-rays, as opposed to time-of-flight delays where the 
higher energy X-rays precede the lower energy X-rays. 

Aschwanden et al. (1997) filtered variations on time scales ~1 s or less out of CGRO 
BATSE flare HXR light curves. They found time delays in the remaining gradually vary- 
ing component to be consistent with magnetic trapping and collisional precipitation of the 
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particles. Trap plasma densities ~ 10*^ cm“^ were deduced. No evidence was found for a 
discontinuity in the delay time as a function of energy and, therefore, for second-step (two- 
stage) acceleration of electrons at energies < 200 keV. 

8.3 Thermal Delays 

The third type, the thermal delay, can be caused by the temperature dependence of cool- 
ing processes, such as by thermal conduction, rdT) a (e.g., Antiochos and Stur- 

rock 1978; Culhane et al. 1994), or by radiative cooling, Xr(T) oc (e.g., Fisher and 
Hawley 1990; Cargill et al. 1995). The observed physical parameters suggest that ther- 
mal conduction dominates in flare loops at high temperatures as observed in soft X- 
ray wavelengths, while radiative cooling dominates in the later phase in postflare loops 
as observed in EUV wavelengths (Antiochos and Sturrock 1978; Culhane et al. 1994; 
Aschwanden and Alexander 2001). When the temperature drops in the decay phase of flares, 
the heating rate can justifiably be neglected and the conductive or radiative cooling rate 
dominate the temperature evolution. Before RHESSI, the cooling curve T (t) in flare plas- 
mas had been studied in only a few flares (e.g., McTiernan et al. 1993; Culhane et al. 1994; 
Aschwanden and Alexander 2001). 

The high spectral resolution of RHESSI data is particularly suitable for any type of ther- 
mal modeling, because we can probe the thermal plasma from ~3 keV up to ~30 keV with 
a FWHM resolution of ~1 keV thanks to the cooled germanium detectors (Lin et al. 2002; 
Smith et al. 2002). This allows us to measure flare temperatures with more confidence. 
A statistical study of flare temperatures measured in the range of J ~ 7-20 MK indeed 
demonstrates some agreement between the values obtained from spectral fitting of RHESSI 
data with those obtained from GOES flux ratios (Battaglia et al. 2005), although RHESSI 
has a bias for the high-temperature tail of the differential emission measure (DEM) distri- 
bution (Aschwanden et al. 2008; Vaananen and Pohjolainen 2007). Of course, we expect an 
agreement between the deduced emission-measure-weighted temperatures only when both 
instruments are sensitive to a temperature range that covers the flare DEM peak. 

A close relationship between the non-thermal and thermal time profiles was found early 
on, in the sense that the thermal emission often closely resembles the integral of the non- 
thermal emission, a relationship that is now known as the Neupert effect (Neupert 1968; 
Hudson 1991; Dennis and Zarro 1993). This relationship is, however, strictly only expected 
for the asymptotic limit of very long cooling times, while a physically more accurate model 
would quantify this effect by a convolution of the non-thermal heating with a finite cool- 
ing time. The deconvolution of the e-folding cooling time in such a model has never been 
attempted statistically and as a function of energy or temperature. Theoretical discussions 
of the Neupert effect, including multiple energy release events, chromospheric evaporation, 
and cooling, can be found in Warren and Antiochos (2004), Liu et al. (2010), and Reeves 
and Moats (2010). 

The cooling time at a given energy can be estimated from the decay time of a flare time 
profile. For instance, the decay times measured with GOES in soft X-rays were found to have 
a median of r^ecay ~ 6 min (Veronig et al. 2002a, 2002b). The observed cooling times have 
typically been found to be much longer than predicted from classical conduction, but shorter 
than the radiative cooling time (e.g., McTiernan et al. 1993; Jiang et al. 2006; Raymond et 
al. 2007). This discrepancy could result from either continuous heating or suppression of 
conduction during the decay phase, or a combination of both. 

The Neupert effect was tested by correlating the soft X-ray peak flux with the (time- 
integrated) hard X-ray flux. A high correlation and time coincidence between the soft X-ray 
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peak and hard X-ray end time was generally found, but a significant fraction of events also 
had a different timing (Veronig et al. 2002c). A delay of 12 s was found in the soft X-ray flux 
time derivative with respect to the hard X-ray flux in SOL2003-1 1-13X05:01 (Ml. 6) (Liu et 
al. 2006, also see Sect. 7.2). Time delays such as this could be related to the hydrodynamic 
flow time during chromospheric evaporation. Tests of the “theoretical Neupert effect,” i.e., 
comparisons of the beam power supply of hard X-ray-emitting electrons and the thermal 
energy of evaporated plasma observed in soft X-rays, found it to strongly depend on the low- 
energy cutoff to the non-thermal electron distribution (Veronig et al. 2005). This provides 
another approach to deducing the energy at which the low-energy cutoff in the electron 
distribution occurs in individual flares. The Neupert effect has also been studied in several 
flares by Ning (2008, 2009), who finds a high correlation between the hard X-ray flux and 
the time derivative of the thermal energy deduced from X-ray spectral fits (Ning 2008) and 
an anti-correlation between the hard X-ray spectral index and the time rate of change of the 
UV flare area observed by TRACE (Ning 2009). 

8.4 Multi-thermal Delay Modeling with RHESSl 


Since major solar flares generally produce a large number of individual postflare loops, giv- 
ing the familiar appearance of loop arcades lined up along the flare ribbons, it is unavoidable 
that each loop is heated up and cools down at different times, so that a spatially integrated 
spectrum always contains a multi-thermal differential emission measure distribution (cf. 
Warren 2006). The resulting multi-thermal bremsstrahlung spectrum (for photon energies e) 
observed in soft X-rays (neglecting the Gaunt factor of order unity). 


/(£) = /o 


/ 


exp(-e/ytflr) dEM(T) 
7W2 df 


dT, 


(8.1) 


is then a function of a multi-thermal differential emission measure distribution dEM(T) = 
n^{T)dV . An example of a multi-thermal spectrum from a differential emission measure 
proportional to up to a maximum temperature of 50 MK is shown in Fig. 14. 

As discussed above, the initial cooling of the hot flare plasma (say at T > 10 MK) is 
generally dominated by conductive cooling (rather than by radiative cooling, which can 
dominate later after the plasma cools to EUV-emitting temperatures of T <2 MK). The 
thermal conduction time has the following temperature dependence: 
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see Aschwanden (2007) for parameter definitions. Since the thermal bremsstrahlung at de- 
creasing photon energies is dominated by radiation from lower temperature flare plasma, the 
conductive cooling time is expected to become longer at lower temperatures (tcond 
Thus, the soft X-ray peak is always delayed with respect to the harder X-ray peaks, reflecting 
the conductive cooling of the flare loops. 

Aschwanden (2007) has measured and modeled this conductive cooling delay Vo„d(e) for 
a comprehensive set of short-duration (< 10 min) flares observed by RHESSl. One example 
is shown in Fig. 15. He finds that the cooling delay At expressed as a function of the photon 
energy e and photon spectral index y can be approximated by 
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Fig. 14 Example of a multi-thermal spectrum with contributions from plasmas with temperatures of 
T = 15, 20 , ... ,50 MK and a DEM distribution of dEM{T)/dT oc T~“^. The individual thermal spectra 
and their sum are shown with thin linestyle, where the sum represents the observed spectrum. Note that the 
photons in the energy range 8 = 5.8-19.4 keV are dominated by temperatures of T = 15-50 MK, which have 
a corresponding thermal energy that is about a factor of (4 -|- 1 /2) = 4.5 lower than the corresponding photon 
energy = 1.3— 4.3 keV). The summed photon spectrum without the high-temperature cutoff approaches 
the power-law function 1(e) a e~^'^ (dotted line) (from Aschwanden 2007). Reproduced by permission of 
the AAS 


(where Xg is the Gaussian width of the time profile peak) and yields a new diagnostic of 
the process of conductive cooling in multi-thermal flare plasmas. In a statistical study of 
65 flares (Aschwanden 2007), 44 (68%) were well fit by the multi-thermal model with 
a best fit value for the exponent of /J = 2.7 ± 1.2, which is consistent with the theo- 
retically expected value of /J = 2.5 according to (8.2). The conductive cooling time at 
To = 11.6 MK (£o = 1 keV) was found to range from 2 to 750 s, with a mean value of 
TcO = 40 s. 

We note that these timing data, as well as thick-target fits to the non-thermal part of spec- 
tra that reveal the evolution of the energy content in accelerated electrons, provide additional 
constraints on models such as the multithread flare model of Warren (2006). 
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Fig. 15 X-ray light curves are shown for SOL2002-02-26T10:27 (C9.6), for energies of 10 keV to 30 keV 
in intervals of 1 keV, observed with RHESSI {left panels). The spectrum is decomposed into thermal and 
non-thermal components (top right panel) and the delay of the peaks at different energies is fitted with a 
thermal conduction cooling time model that has a scaling of ^ (right bottom panel). The 

best fit shows a power index of ^ = 2.8, which is close to the theoretically expected value of ^ = 5/2 
(8.2). The full delay of the thermal component is indicated with a thin curve (bottom left panel), while the 
weighted (thermal-|-non-thermal) fit is indicated with a thick curve (from Aschwanden 2007). Reproduced 
by permission of the AAS 
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9 Hard X-ray Spectral Evolution in Flares 

9.1 Observations of Spectral Evolution 

The non-thermal hard X-ray emission from solar flares, best observed in the 20 to 100 keV 
range, is highly variable. Often several emission spikes with durations ranging from seconds 
to minutes are observed. In larger events, sometimes a more slowly variable, long duration 
emission can be observed in the later phase of the flare. Hence, most flares start out with an 
impulsive phase, while some events, mostly large ones, show the presence of a late gradual 
phase in the hard X-ray time profile. 

While these two different behaviors can already be spotted by looking at lightcurves, 
they also are distinct in their spectral evolution. The impulsive spikes tend to be harder at 
the peak time, and softer both in the rise and decay phase. The spectrum starts soft, gets 
harder as the flux rises and softens again after the maximum of the emission. This pattern of 
the spectral evolution is thus called soft-hard-soft (SHS). On the other hand, in the gradual 
phase, the flux often slowly decreases, while the spectrum stays hard or gets even harder. 
This different kind of spectral evolution is called soft-hard-harder (SHH). 

Historically, both the SHS (Parks and Winckler 1969; Kane and Anderson 1970) and 
the SHH behavior (Frost and Dennis 1971) were observed in the early era of hard X-ray 
observations of the Sun. Subsequent investigation confirmed both the SHS (Benz 1977; 
Brown and Loran 1985; Lin and Schwartz 1987; Gan 1998; Fletcher and Hudson 2002; 
Hudson and Farnlk 2002) and the SHH (Cliver et al. 1986; Kiplinger 1995; Saldanha et al. 
2008; Grigis and Benz 2008) behavior. The SHH behavior has been found to be correlated 
with proton events in interplanetary space (Kiplinger 1995; Saldanha et al. 2008; Grayson 
et al. 2009). 

Evidence for hard-soft-hard (HSH) spectral evolution at energies above ~50 keV has 
been reported for multiple spikes in SOL2004-1 1-03T03:35 (Ml. 6) (Shao and Huang 
2009b). SHS behavior was observed at lower energies. This HSH behavior might be ex- 
plained by albedo, which typically peaks around 30-40 keV (see Kontar et al. 2011), but the 
authors corrected for albedo from isotropically emitted photons. A likely explanation is that 
the spikes overlie a harder, gradually varying component, possibly emission from trapped 
electrons (Sect. 8.2). 

While all these observations established the qualitative properties of the spectral evolu- 
tion, a statistical analysis of the quantitative relation between the flux and spectral index had 
not been performed in the pre-RHESSI era. Here, we summarize RHESSI results investigat- 
ing quantitatively the spectral evolution of the non-thermal component of the hard X-ray 
emission, as well as the theoretical implications. More details can be found in Grigis and 
Benz (2004, 2005, 2006). 

To quantify the spectral evolution, a simple parameterization for the shape of the non- 
thermal spectrum is needed. Luckily, in solar flares the spectrum is well described by a 
power law in energy, which often steepens above 50 keV. Such a softening of the spectrum 
can be modeled by a broken power-law model. However, it is difficult to observe such a 
downward bending at times of weak flux, because the high-energy region of the spectrum is 
lost in the background. As a compromise, Grigis and Benz (2004) fitted the data to a single 
power-law function at all times. Although the single power law does not always provide 
a good fit to the spectra, it provides a characteristic spectral slope and ensures a uniform 
treatment of the spectra at different times. 

The two free parameters of the power-law model are the spectral index y and the power- 
law normalization at the reference energy Sq. The reference energy eo is arbitrary, but 
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Fig. 16 Time evolution of the spectral index y {upper curve, linear scale on right) and the flux normaliza- 
tion 735 {lower curve, logarithmic scale on left) of the non-thermal component in SOL2002-1 1-09X13:23 
(M4.9). Different emission spikes are shown in different colors (after Grigis and Benz 2004). Reproduced 
with permission ©ESO 


fixed, usually near the logarithmic mean of the covered energy range. In the RHESSI spectral 
analysis software, OSPEX, Sq = 50 keV by default. The time dependent spectrum is given 
by 

/eV’'” 

/(£,0 = 4o(0(-j • (9.1) 

A representative sample of 24 solar flares of GOES magnitudes between Ml and XI 
was selected by Grigis and Benz (2004). The spectral model (9.1), with the addition of an 
isothermal emission component at low energies, was fitted with a cadence of one RHESSI 
spin period (about 4 s). This delivered a sequence of measurements of the quantities /e„(f) 
and y (f) for each of the 24 events, covering a total time of about 62 minutes of non-thermal 
hard X-ray emission. For these events, Eq = 35 keV was chosen, a meaningful energy which 
lies about in the middle of the range where the non-thermal emission is best observed in 
these M-class flares. 

An example of the measured time evolution of the spectral index y and the flux normal- 
ization 735 for the longer-lasting event of the set is shown in Fig. 16. A correlation in time 
between the two curves can be readily seen. Single emission spikes are plotted in differ- 
ent colors, so that the soft-hard-soft evolution can be observed during each spike (with the 
exception of the late, more gradual phase, where the emission stays hard as the flux decays). 

As there is an anti-correlation in time between log hs (f) and / (f), a plot of one parameter 
as a function of the other, eliminating the time dependence, shows the relationship between 
them. Figure 17 shows plots of y vs. h5 for 3 events where there are only one or two 
emission peaks. The points in the longer uninterrupted rise or decay phase during each event 
are marked by plus symbols. A linear relationship between log hs and y can be seen during 
each phase, although it can be different during rise and decay. 

On the other hand, a plot of all the 911 fitted model parameters for all the events show 
a large scatter, as shown in Fig. 18. The large scatter can be understood as originating from 
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Fig. 17 Spectral index y vs. flux normalization 735 for three events, showing the linear dependence of single rise and decay phases of emission spikes on a log-linear scale. 
Dots mark results from individual spikes, while pluses mark the longer rise or decay phase (from Grigis and Benz 2004). Reproduced with permission ©ESO 
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Fig. 18 Plot of the spectral 
index y versus the fitted 
non-thermal flux at 35 keV (given 
in photons s“* cm“^ keV“*). 

All 911 data points from the 
24 events are shown (from Grigis 
and Benz 2004). Reproduced 
with permission ©ESO 
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the superposition of data from a large numbers of different emission spikes, each featuring 
linear trends with different parameters. This plot does demonstrate, however, the tendency 
for flatter spectra to be associated with more intense flares. 

RHESSI observations of the gradual phase of large solar flares (Grigis and Benz 2008) 
and its relation with proton events (Saldanha et al. 2008; Grayson et al. 2009) have shown 
that the hardening behavior is complex and cannot be characterized by a continuously in- 
creasing hardness during the event. Therefore the soft-hard-harder (SHH) denomination 
does not accurately reflect the observed spectral evolution. Rather, phases of hardening (or 
even approximatively constant hardness) are often seen in larger events as the flux decays 
(Kiplinger 1995). The start of the hardening phase can happen near the main peak of the 
flare, or later. The end of hardening can even be followed by new impulsive SHS peaks. The 
most recent statistical study of the correlation of SHH behavior with proton events (Grayson 
et al. 2009) found that in a sample of 37 flares that were magnetically well-connected to 
Earth, 18 showed SHH behavior and 12 of these produced solar energetic particle (SEP) 
events. None of the remaining 19 flares that did not show SHH behavior produced SEP 
events. 

9.2 Interpretation of Spectral Evolution 

Can we explain the soft-hard-soft spectral behavior theoretically? The problem here is that 
many effects contribute to the properties of the high-energy electron distribution whose 
bremsstrahlung hard X-rays are observed by RHESSI and similar instruments. We can iden- 
tify three main, closely related classes of physical processes that affect the distribution of 
the electrons and the spectrum of the X-ray photons they generate: (1) the acceleration of 
part of the thermal ambient plasma, (2) the escape from the acceleration region, and (3) the 


^ Springer 


G.D. Holman et al. 


transport to the emitting region. The photon spectrum also depends on the properties of the 
hremsstrahlung emission mechanism. 

Miller et al. (1996) proposed a stochastic acceleration mechanism where electrons are 
energized by small- amplitude turbulent fast-mode waves, called the transit-time damping 
model. They showed that their model could successfully account for the observed number 
and energy of electrons accelerated above 20 keV in subsecond spikes or energy release 
fragments in impulsive solar flares. However, they made no attempt to explain the observed 
hard X-ray spectra (which are softer than predicted by the transit-time damping model) and 
did not consider spectral evolution. Furthermore, this approach does not account for particle 
escape. Grigis and Benz (2006) extended the model with the addition of a term describing 
the escape of the particles from the acceleration region, as in the model of Petrosian and 
Donaghy (1999). To ensure conservation of particles, they also add a source term of cold 
particles coming into the accelerator (such as can be provided by a return current). 

The stochastic nature of this acceleration model implies that the electrons undergo a 
diffusion process in energy space. Mathematically, the acceleration is described by the fol- 
lowing convective-diffusive equation: 

? = T l^t^'OcoLL + Dj)f] - ^[(Acoll + At)/] - S(E) • / + Q{E), (9.2) 

dt 2 dE^ dE 


where f(E) is the electron density distribution function, Dj and Aj are, respectively, the 
diffusion and convection coefficients due to the interactions of the electrons with the acceler- 
ating turbulent waves, Dcoll and Acoll are, respectively, the diffusion and convection coef- 
ficients due to collisions with the ambient plasma, S{E) is the sink (escape) term, and Q(E) 
is the source (return current) term. The escape term is proportional to v{E) /r, where v{E) is 
the electron speed, and r is the escape time. The escape time can be energy-dependent, but 
for simplicity it is initially kept constant. The longer the escape time, the better the particles 
are trapped in the accelerator. The source term is in the form of a Maxwellian distribution 
of electrons with the same temperature as the ambient plasma. 

The coefficients Dj and Aj are proportional to the dimensionless acceleration parameter 
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where Uj and are, respectively, the energy densities of the turbulent waves and of the 
ambient magnetic field, (k) is the average wave vector, and is the proton gyrofrequency. 

Equation (9.2) can be solved numerically until an equilibrium state (df/dt = 0) is 
reached. The equilibrium electron spectra from the model are controlled by two parameters: 
the acceleration parameter 7 acc described above and the escape time r. Above 10-20 keV, 
the collision and source terms in (9.2) can be neglected, since they apply to the ambient 
Maxwellian, and thus the equilibrium spectra depend to a first approximation only on the 
product / = I ACC ■ T- 

Figure 19 shows the equilibrium electron spectra for different values of / = Iacc ■ t. As 
/ increases, the spectrum gets harder and harder. To explain the soft-hard-soft effect, either 
the acceleration or the trapping efficiency (or both) must increase until the peak time, and 
then decrease again. We note that this model does not include magnetic trapping (other than 
in the magnetic turbulence itself), which can alter the computed electron spectra and their 
time evolution (e.g., Metcalf and Alexander 1999). 

To see whether this produces the linear relation between the spectral index and the log 
of the flux normalization, Grigis and Benz (2006) computed the hard X-ray emission from 
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Fig. 19 Accelerated electron density distributions with different values of the power-law index resulting 
from changes in = I\cc ' ^ ■ The dashed curve represents the ambient Maxwellian distribution. The two 
dolled lines indicate the energy range used for the computation of the power-law index S shown above each 
spectrum. Harder spectra have a larger /j value (from Grigis and Benz 2006). Reproduced with permission 
©ESO 

these model electron spectra. Since these are equilibrium spectra, thin-target emission was 
computed. They then plotted the spectral index vs. the flux normalization of the resulting 
photon spectra. Since the spectra are not power-law, but bend down, they fit a power-law 
model to the model photon spectrum in a similar range as the one used for the observations. 

Figure 20 shows the computed values for the spectral indices and flux normalizations for 
both the electron and the photon spectrum from the model. The results show that there is 
indeed a linear relation between the spectral index and the log of the flux normalization. 

An alternative mechanism that could be responsible for soft-hard-soft spectral evolution 
is return current losses as the electrons propagate to and within the thick-target footpoints 
of the flare loop (Zharkova and Gordovskyy 2006). The highest electron energy to which 
return current losses are significant is proportional to the return current electric field strength, 
which is in turn proportional to the electron beam flux density (see Sect. 5). Therefore, 
as the electron flux density increases and then decreases, the low-energy part of the X-ray 
spectrum flattens to higher and then lower energies as the return current electric field strength 
increases and then decreases. The net effect is SHS spectral evolution below the maximum 
energy for which return current losses are significant during the flare. The observation of 
SHS behavior in coronal X-ray sources, however, indicates that this spectral evolution is a 
property of the acceleration mechanism rather than a consequence of energy losses during 
electron propagation (Battaglia and Benz 2006, see Sect. 10.3). 
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Fig. 20 Model results for the spectral index and flux normalization for electrons and photons. The dashed 
line is the best straight-line fit to the model results (in the range of spectral indices from 2 to 8 for the 
electrons, and 3 to 9 for the photons), corresponding to a pivot-point behavior (from Grigis and Benz 2006). 
Reproduced with permission ©ESO 


Are there two stages of electron acceleration, one responsible for the impulsive phase 
and one for the gradual phase? RHESSI spectroscopy and imaging of a set of 5 flares with 
hardening phases (Grigis and Benz 2008) showed that there is no discontinuity in the motion 
of footpoints at the onset of hardening and no clear separation between the impulsive and 
the gradual phase: the former seems to smoothly merge into the latter. This supports the 
view that the same acceleration mechanism changes gradually in the later phase of the flare, 
rather than a two stage acceleration theory. The hardening phase may in fact be caused by 
an increase in the efficiency of trapping of the electrons above 100 keV. 

The underlying cause of the SHS spectral evolution has been addressed in terms of the 
stochastic acceleration model by Bykov and Fleishman (2009) and Liu and Fletcher (2009). 
Bykov and Fleishman consider acceleration in strong, long-wavelength MFID turbulence, 
taking into account the effect of the accelerated particles on the turbulence. They argue 
that the electron spectrum flattens during the linear acceleration phase, while the spectrum 
steepens during the nonlinear phase when damping of the turbulence because of the parti- 
cle acceleration is important, giving SHS spectral evolution. They argue that SHH evolu- 
tion will be observed when the injection of particles into the acceleration region is strong. 
Liu & Fletcher also argue that the SHS evolution results from dependence of the electron 
distribution power-law index on the level of turbulence as it increases and subsequently 
decreases. They attribute changes in the SHS correlation during a flare to changes in the 
background plasma, likely due to chromospheric evaporation. 

We note that simple direct-current (DC) electric field acceleration of electrons out of the 
thermal plasma can produce the SHS spectral evolution. The flux of accelerated electrons 
and the maximum energy to which electrons are accelerated and, therefore, the high-energy 
cutoff to the electron distribution, increase and decrease together as the electric field strength 
increases and decreases (Holman 1985). The X-ray spectrum is steeper at energies within 
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one to two orders of magnitude below the high-energy cutoff (Holman 2003). In large flares, 
however, where the X-ray spectrum is observed to continue to MeV energies or higher, there 
is no evidence for a high-energy cutoff in the appropriate energy range. Therefore, at least for 
large flares with spectra extending to high energies, a simple DC electric held acceleration 
model does not appear to be appropriate. 


10 The Connection Between Footpoint and Coronal Hard X-ray Sources 

Hard X-ray (HXR) sources at both footpoints of a coronal loop structure have been observed 
since Hoyng et al. (1981). As reviewed in Sects. 1 and 2, they are understood to be thick- 
target bremsstrahlung emission produced by precipitating electrons, accelerated somewhere 
in or above the loop. A third HXR source situated above the looptop (see Krucker et al. 
2008a, for a review) was first noted by Masuda et al. (1994) in Yohkoh observations. The 
nature of this coronal HXR source has remained uncertain, but in simple solar flare models 
with reconnection and particle acceleration in the corona, we expect some relation between 
coronal HXR sources and footpoints. RHESSI has enabled us to study events featuring coro- 
nal HXR sources and footpoints simultaneously. By studying the behavior of the sources in 
time and the relations between them, we can address questions such as: Are both coronal and 
footpoint emissions caused by the same electron population? How is such an electron beam 
modified in the loop (collisions, return currents, trapping, etc.)? Is SHS behavior (Sect. 9.1) 
a transport effect produced by collisions or return currents, or is it a feature imposed by the 
acceleration mechanism? 

10.1 Imaging Spectroscopy 

RHESSI has provided the possibility of obtaining simultaneous, high-resolution imaged 
spectra at different locations on the Sun. One can therefore study each source separately 
in events with several contemporaneous HXR sources. The high spectral resolution has al- 
lowed a reliable differentiation between thermal and non-thermal emission to be made in 
many flares. Furthermore, RHESSEs imaging spectroscopy has allowed differences in indi- 
vidual flare source spectra and their evolution to be studied in considerable detail. 

Imaged spectra and the relative timing of sources in three flares, including the limb flare 
SOL2002-02-20T1 1 :07, were studied by Krucker and Lin (2002). Sui et al. (2002) analyzed 
and modeled the two footpoint sources and a high, above-the-looptop hard X-ray source 
observed in this flare. Emslie et al. (2003) analyzed SOL2002-07-23T00:35 (X4.8) flare 
with four HXR sources observed by RHESSI. They found a coronal source with a strong 
thermal component, but the non-thermal component could not be studied due to severe pulse 
pile-up. Battaglia and Benz (2006) studied five M-class events. Due to the smaller pile-up 
amount in those events, studying the non-thermal coronal emission was possible. The results 
of these studies are summarized below. 

10.2 Relation Between Coronal and Footpoint Sources 

The quantitative relations between the footpoints and the coronal source and between the 
two footpoints can give information about the physical mechanisms at work in a solar flare. 
Simple models envision a beam of accelerated electrons encountering a low-density region 
in the corona, leading to thin-target bremsstrahlung. When the same electron beam reaches 
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the chromosphere, the particles are fully stopped in the dense material, producing thick- 
target emission. Assuming an electron power-law distribution for the electron energy E of 
the form 

^(E) = AE~^ ( 10 . 1 ) 

producing thin-target bremsstrahlung in the coronal source, the observed photon spectrum 
has spectral index = S -|- 1 (2.10). Reaching the chromosphere, the accelerated elec- 
trons will be fully stopped, producing thick-target bremsstrahlung with a photon spectral 
index y,hick = 3—1 (2.12). In such a simple scenario one would therefore expect a differ- 
ence in the photon spectral index y,/„„ — y,hick = 2 between the coronal source and the foot- 
points. Further, the two footpoints should be of equal hardness and intensity if one assumes 
a symmetric loop and symmetric injection of particles into the legs of the loop. 

10.2.1 Observed Difference Between Coronal and Eootpoint Spectral Indices 

A sample of flares observed with Yohkoh to have coronal HXR sources was studied by 
Petrosian et al. (2002). They found that the spectral index of the coronal sources was, on 
the average, steeper by 1 than the spectral indices of the footpoint sources. Sui et al. (2002) 
also found a spectral index difference of 1 for SOL2002-02-20T1 1:07 (C7.5) observed with 
RHESSl. 

Battaglia and Benz (2006) found that the coronal source was softer than both footpoints 
for all of their five events in nearly all analyzed time bins. Figure 21 (top left) shows an 
image of SOL2005-07-13T14:49 (M5.0) in the 34-38 keV energy band. The two footpoints 
are visible, as well as the 50 and 80% contours of the coronal source taken from a 10- 
12 keV image. Spectra and spectral fits are shown for the two footpoints and the coronal 
source. The steepness of the coronal source spectrum (number 3 in the figure) relative to 
the spectra from the footpoints is apparent. However, the quantitative difference between 
the values of the spectral index obtained for the coronal source and the footpoints often 
differed significantly from 2. For the five Hares analyzed, the smallest mean difference in 
the spectral indices, averaged over time, was 0.59±0.24. The maximum mean difference, 
averaged over time, was 3.68±0.14. These clearly contradict the theoretical expectation 
summarized above. Simple thin-thick target scenarios do not seem to work in most cases 
and additional effects need to be considered. 

Evidence for two populations of coronal source non-thermal spectra was found by Shao 
and Huang (2009a). They compare coronal and footpoint spectral indices at 28 hard X-ray 
peaks from 13 single-loop flares observed by RHESSl. The spectral index in the coronal 
sources was determined from an isothermal plus power-law fit below 30 keV, while the 
footpoint spectral indices were determined from a power-law fit at 30-60 keV photon ener- 
gies. They argue that the coronal spectra can be divided into two groups. One, for which the 
coronal spectral index is greater than 5, is well correlated with the footpoint spectral index, 
and the difference in the indices ranges from 2^. For the other, where the spectral indices 
are anticorrelated, the coronal spectral index is less than 5, and the difference in the indices 
ranges from 0-2. For the group of anticorrelated spectral indices, the coronal spectral index 
is correlated with the photon flux, while the footpoint spectral index is anticorrelated with 
the photon flux for both groups. These are intriguing results if confirmed by future studies. 

10.2.2 Differences Between Footpoints 

No significant difference was found in the spectral indices for the two footpoints in 
SOL2002-02-20T11:49 (C7.5) by Krucker and Lin (2002) and Sui et al. (2002). Piana et 
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RHESSI 34-38 keV 13-Jul-2005 14:14:00.000 UT 



Fig. 21 Top left. Composite CLEAN image of a RHESSI event with three hard X-ray sources. The footpoints 
(labeled 1 & 2) are visible on the solar disc in an image made at 34—38 keV. The position of the coronal source 
(labeled 3) high above the limb is indicated by the 50 and 80% white contours taken from a 10-12 keV image. 
Plots 1-3 show spectra and normalized residuals over the fitted energy range for the north footpoint (1), south 
footpoint (2), and coronal source (3) (after Battaglia and Benz 2006). Reproduced with permission ©ESO 


al. (2007) inverted count visibility spectra for this flare to obtain mean electron flux distri- 
butions for the footpoints. They found the mean electron flux distribution function at the 
northern footpoint to be somewhat steeper (A3 ~ 0.8) than that derived for the southern 
footpoint. They also found the distribution function for the region between the footpoints 
(not the coronal source studied by Sui et al.) to be steeper than the footpoint distribution 
functions (A3 ~ 1.6 relative to the southern footpoint) and to substantially steepen at ener- 
gies above ~60 keV. 

Krucker and Lin (2002) found that, when a connection between footpoints could be de- 
termined, the footpoints brightened simultaneously (to within the ~1 s time resolution of 
the observations) and had similar spectra. 

Differences of 0. 3-0.4 between the spectral indices of two footpoints in SOL2002-07- 
23T00:35 (X4.8) were reported by Emslie et al. (2003). 

For the flares analyzed by Battaglia and Benz (2006), a significant difference was found 
in only one out of five events. For all other events, the mean difference in yfp was zero within 
the statistical uncertainty. 
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Different spectra at the two footpoints imply an asymmetric loop. Such an asymmetry 
can result, for example, from different column densities or different beam fluxes and corre- 
sponding return current energy losses in the legs of the loop. It could also result from asym- 
metric magnetic trapping within the loop (e.g., Alexander and Metcalf 2002). In a study 
of 53 flares showing two HXR footpoints, Saint-Hilaire et al. (2008) found that footpoint 
asymmetry was greatest at the time of peak HXR flux and the difference in the footpoint 
spectral indices Ay rarely exceeded 0.6. In most cases they found the footpoint asymmetry 
to be inconsistent with different column densities in the two legs of the loops. 

In SOL2003- 10-29X20:49 (XIO.O) Liu et al. (2009a) found that the brighter HXR 
footpoint was marginally, but consistently harder than the dimmer footpoint by A/ = 
0.15 ± 0.13. They concluded that neither asymmetric magnetic mirroring nor asymmet- 
ric column density alone can explain the full time evolution of the footpoint HXR fluxes 
and spectral indices. However, a self-consistent explanation might be obtained by consider- 
ing these two effects together and/or in combination with one or more additional transport 
effects, such as nonuniform target ionization, relativistic beaming, and return current losses. 

10.3 Spectral Evolution in Coronal Sources 

Previous observations of SHS spectral evolution (see Sect. 9.1) were made with full-Sun 
spectra which, except for over-the-limb events, are typically dominated by footpoint emis- 
sion. Battaglia and Benz (2006), in their imaging spectroscopy study, found that the coronal 
source itself shows SHS evolution. This is illustrated in Fig. 22. This finding implies that 
SHS is not caused by transport effects within the flare loop, but is rather a property of 
the acceleration mechanism itself. Indeed, Grigis and Benz (2006) showed that SHS can 
be reproduced for electron spectra in a transit-time-damping, stochastic-acceleration model 
(Sect. 9.2). 

10.4 Interpretation of the Connection Between Footpoints and The Coronal Source 

In the above account, emphasis was given to the difference in the spectral index between the 
coronal source and footpoints. Assuming a thin target in the corona and a thick target at the 
footpoints, one would expect a difference of two. However, whether the coronal source acts 
as thin- or thick-target depends on the energy of the accelerated electrons and the column 
density in the corona. Veronig and Brown (2004), for example, found coronal sources with 
column densities high enough to act as thick targets for electrons with energies up to 60 keV. 

As early as 1976, Melrose and Brown (1976) showed that magnetic trapping with colli- 
sional scattering of electrons out of the trap can lead to a thick-target coronal source. The 
coronal source transitions through a thin-thick period, with the time scale for this transition 
depending on the electron energy and the plasma density in the trap. The trapping essen- 
tially increases the effective column density in the corona. Metcalf and Alexander (1999) 
analyzed six flares with coronal sources observed by Yohkoh and found that three of the six 
flares showed properties consistent with trapping. 

A simple 1-D model that described the coronal emission as intermediate thin-thick, 
depending on electron energy, was developed by Wheatland and Melrose (1995). In this 
model a high-density region (>10*^ cm“^) is hypothesized to be present at or above 
the top of the flare loop. The model makes predictions for the shape of the coronal and 
footpoint spectra and the relations between them. Fletcher (1995) obtained Monte Carlo 
solutions to the Fokker-Planck equation to show that, with the inclusion of high elec- 
tron pitch angles and collisional scattering, a compact coronal X-ray source is produced 
at the top of a loop with a constant coronal density ~3 x 10*° cm“^. Holman (1996) 
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Fig. 22 Top: GOES 1-8 A light curve of SOL2003- 10-24X02:54 (M7.6). Middle: RHESSI 25-50 and 
50-100 keV light curves near the peak of the GOES flare. Bottom: time evolution of fitted coronal source 
flux at 35 keV (F 35 , * symbols, left log scale) and spectral index (y, + symbols, right log scale) displaying 
SHS evolution (after Battaglia and Benz 2006). Reproduced with permission ©ESO 


showed that, even in the simple 1-D model, a compact coronal source is produced when 
electrons are injected into a loop with a constant coronal density '^2 x 10^^ cm“^ (see 
hesperia.gsfc.nasa.gov/sftheory/loop.htm). A compact coronal HXR source can also be pro- 
duced if there is a compact magnetic trap at or above the top of the loop. Fletcher and 
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Fig. 23 Left: spectra for coronal source {red) and footpoints (blue) according to the model of Wheatland and 
Melrose (1995). The spatially integrated spectrum is shown in violet. Right: observed RHESSI spectra for the 
event SOL2003-10-24T02:54 (M7.6). Isothermal and power-law hts to the coronal (crosses) and footpoint 
(dots) spectra are shown. The vertical line indicates the predicted critical energy for the transition between 
thin and thick target (after Battaglia and Benz 2007). Reproduced with permission ©ESO 


Martens (1998) showed that, with such a trap, a significant coronal X-ray source can be pro- 
duced at plasma densities as low as ~4 x 10^ cm“^. Petrosian and Donaghy (1999) showed 
that the coronal HXR source can be a consequence of acceleration and trapping by turbu- 
lence or plasma waves. In their stochastic acceleration model, the difference between the 
coronal and footpoint spectra is explained by the energy-dependent time scale for electrons 
to escape the acceleration region. 

The left panel of Fig. 23 illustrates the model of Wheatland and Melrose (1995). The 
spatially integrated spectrum (violet) is the power-law spectrum (thick-target, Ythick = 3 — 1) 
expected for a single-power-law electron distribution with no low- or high-energy cutoffs 
and no thermal component. For e = y/2KN (see (2.4)), the spectrum is dominated 
by thick- target radiation from the coronal source (red). There is a low-energy cutoff in the 
electron distribution at the footpoints at E ^ ^2KN because of the energy losses in the 
coronal source. The spectrum is dominated by thick-target radiation from the footpoints 
(blue) where e )^> £<,. It is in this regime that the radiation from the coronal source is thin- 
target and the spectral index of the coronal source is steeper by 2 than that of the footpoints. 
These spectra are characteristic of all the models reviewed above. 

Sui et al. (2002) compared the RHESSI observations of SOL2002-02-20T1 1:07 (C7.5) 
to a model with a constant-coronal-density loop and no magnetic trapping. They used a 
finite difference method (e.g., McTiernan and Petrosian 1990; Holman et al. 2002) to ob- 
tain steady-state solutions to the Fokker-Planck equation with collisional scattering and en- 
ergy losses. Model images were convolved with the RHESSI response to produce simulated 
RHESSI observations for direct comparison with the SOL2002-02-20T1 1:07 flare images 
and imaged spectra. They found that, after obtaining a power-law model spectrum with an 
index of / = 3 that agreed with the observed footpoint spectra, the effective spectral index 
of the coronal source from the model (/ = 4.7) was significantly steeper than that obtained 
for the flare (y =4). 

Battaglia and Benz (2007) compared the model of Wheatland and Melrose (1995) to 
the results of their study of five flares observed by RHESSI. The right panel of Fig. 23 
shows observed spectra and spectral fits for one particular event. The observed spectra were 
dominated by thermal coronal emission at low energies. Therefore, not all of the model 
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Fig. 24 SOL2003-10-28T11:10 (X17.2). Left, bottom: (see Wai'muth et al. 2007) 200-400 MHz radio spec- 
trum (Astrophysical Institute Potsdam) showing the signature of the outflow termination shock (TS, starting at 
1 1 :02:47 UT). Left, top: INTEGRAL count rates at 150 keV and 7.5-10 MeV. Right: (after Aurass et al. 2007, 
reproduced with permission ©ESO.) radio source positions (Nan9ay Radio Heliograph, 327 MHz) overlaid 
on a SOHO-Hn image (11:47 UT 195 A). The bright areas are EUV flare ribbons in ARI0486. RHESSI 
HXR centroids are shown as The integration time intervals are: for the TS source SW of AR10486 
11:02:45-11:03:15 UT, for the continuum source CONT N of AR10486 11:13-11:17 UT, respectively (see 
also Eig. 25). The radio contours are at 50, 70, and 99.5% of the peak flux value 


predictions could be tested. However, the observed relations between the spectra did not 
agree with the predictions of the model. For the flare in Fig. 23, for example, the difference 
between the coronal source and footpoint spectral indices at the higher photon energies 
is 3.8 ±0.1, not 2. Also, an estimate of the column density in the coronal source gives 
^2KN ~ 10-15 keV, while the intersection of the coronal and footpoint spectra is found 
to be at e ~ 23 keV. Battaglia and Benz (2008) have found that this large difference in the 
spectral indices is consistent with spectral hardening caused by return current losses (see 
Sect. 5). 


11 Identification of Electron Acceleration Sites from Radio Observations 

While energetic electrons excite hard X-ray emission during their precipitation into the 
dense layers of the solar atmosphere, they can also excite decimeter and meter wave ra- 
dio emission during propagation and trapping in magnetic field structures in the dilute solar 
corona. The radio emission pattern in dynamic spectrograms can give information about the 
electron acceleration process, the locations of injection of electrons in the corona, and the 
properties of the coronal magnetoplasma structures. 

Here we take as an example SOL2003-10-28T1 1:10 (X17.2) (shown in Fig. 24). Differ- 
ent acceleration sites can be discriminated during the impulsive and the gradual flare phases. 
Radio spectral data from the Astrophysical Institute Potsdam (AIP; Mann et al. 1992), imag- 
ing data from the Nan 9 ay Radio Heliograph (NRH, Kerdraon and Delouis 1997), and hard 
X-ray (RHESSI, INTEGRAL) data were combined in the analysis of this event. The con- 
clusion was reached that a nondrifting, high-frequency type II radio burst signature in the 
radio spectrum coincided with a powerful electron acceleration stage. Simultaneously with 
the nondrifting type II signature, highly relativistic (>10 MeV) electrons were observed in 
the impulsive phase of the flare (Fig. 24, upper left). The radio spectrum suggests that this 
can be due to acceleration at the reconnection outflow termination shock (Aurass and Mann 
2004), as predicted for a classical two-ribbon flare (Forbes 1986; Tsuneta and Naito 1998; 
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Fig. 25 Timing of the source 
CONT in Fig. 24: the NRH 
327 MHz flux curve (asterisks) 
versus the GOES 0.5^ A flux 
curve (solid line, partly 
off-scale). Insef. SOHOfElT 
image showing the radio source 
centroid (white asterisk) and 
RHESSI HXR centroids as in 
Fig. 24. Thick bar. the start time 
of GeV-energy proton injection 
in space (after Aurass et 
al. 2007). Reproduced with 
permission ©ESO 



Aurass et al. 2002). The radio source site is observed about 210 Mm to the SW of the flaring 
active region (TS in Fig. 24, right). In this direction, TRACE and 50//0/LASC0' C2 im- 
ages reveal dynamically evolving magnetoplasma structures in an erupting arcade (Aurass 
et al. 2006). For realistic parameters derived from these observations (the geometry, density, 
temperature, and low magnetic held values of ~5 Gauss), Mann et al. (2006) have found 
that a fully relativistic treatment of acceleration at the fast-mode outflow shock can explain 
the observed fluxes of energetic particles (see Zharkova et al. 201 1). 

In the main flare phase of the same event, an additional radio source (CONT in Fig. 24) 
was found, lasting for ~10 min, indicating the presence of another acceleration site. No 
X-ray, EUV, or Ha emission was observed at the location of this radio source. Figure 25 
gives the timing and the source position with respect to the flaring active region. CONT is a 
m-dm-continuum source with fiber burst fine structure. Fiber bursts are excited by whistler 
waves propagating along held lines of the coronal magnetic held. As marked by a bold bar 
in the figure, the time of the CONT emission is also the start time of GeV proton injection 
in space. Aurass et al. (2006) have shown that this source site is not far from an open field 
(particle escape) region in the potential coronal magnetic field. The source briefly flashes up 
already in the early impulsive phase. Based on a new method of fiber burst analysis (Aurass 
et al. 2005; Rausche et al. 2007), Aurass et al. (2007) argue that this source most likely 
indicates acceleration at a contact between separatrix surfaces of different magnetic flux 
systems. 

Radio observations of flares and their implications are further addressed in White et al. 

(2011). 


12 Discussion and Conclusions 

12. 1 Implications of X-ray Observations for the Collisional Thick- Target Model 

As discussed in Sect. 2, the core assumption of the collisional thick-target model is that 
the spatially integrated hard X-ray emission from non-thermal electrons is bremsstrahlung 

' Large Angle and Spectrometric COronagraph. 
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(free-free radiation) from electrons that lose all their suprathermal energy through colli- 
sional losses in the ambient plasma as they simultaneously radiate the hard X-rays. “Simul- 
taneously” means within the observational integration time. This implies that all electrons 
that contribute significantly to the observed radiation reach a plasma dense enough or, more 
precisely, traverse a high enough column density for all of their suprathermal energy above 
the observed photon energies to be collisionally lost to the ambient plasma within the inte- 
gration time. For typical > 1 s integration times, these conditions are met when the electrons 
stream downward from the corona into the increasingly dense plasma of the solar transition 
region and chromosphere. 

Since the thick-target model is often implicit in our interpretation of the hard X-ray emis- 
sion from flares, it is important to keep the underlying assumptions in mind and test the 
model while at the same time applying it to flare observations. We have discussed above 
several physical processes that, if significant, change the conclusions of the simple colli- 
sional thick-target model regarding the electron distribution produced in the acceleration 
region. These processes occur in either the thick-target region itself, or during the propa- 
gation of the electrons from the acceleration region to the thick-target region. Only with 
the high spectral resolution and imaging of RHESSI has it become possible to observation- 
ally address these processes. Even with the RHESSI observations, however, it is difficult to 
conclusively determine the importance of each process. 

A physical process that distorts the emitted X-ray spectrum is albedo (Sect. 3.4 and Kon- 
tar et al. 2011). Fortunately, the albedo contribution to the X-ray spectrum can be easily cor- 
rected on the assumption that the X-ray photons are isotropically emitted. This correction is 
available in the RHESSI spectral analysis software. If the photons are significantly beamed 
downward, however, the distortion of the spectrum can be substantially greater than that 
from isotropically emitted photons. An anisotropic photon distribution results from emit- 
ting electrons with an anisotropic pitch-angle distribution. The degree of anisotropy of the 
electron pitch-angle distribution also quantitatively affects conclusions from the thick-target 
model concerning the acceleration process. Therefore, it is important to better determine the 
pitch-angle distribution of the emitting electrons and the contribution of albedo to the hard 
X-ray spectrum (see Kontar et al. 2011). 

The simple collisional thick-target model assumes that the target plasma is fully ionized. 
We have seen, however, that a nonuniformly ionized target region can produce an upward 
kink, or “chicane,” in an otherwise power-law X-ray spectrum (Sect. 4). This spectral shift 
can provide a valuable diagnostic of the ionization state of the target plasma and its evo- 
lution. It is likely, however, that the power-law spectrum below the chicane is hidden by 
thermal radiation. The chicane is then observed only as a downward break in the spectrum 
at energies above those dominated by the thermal emission. The upper limit on the magni- 
tude of the break provides a method for ruling out nonuniform ionization as the sole cause 
of large spectral breaks. To further distinguish this break from spectral breaks with other 
causes, it is important to better determine the degree of ionization as a function of column 
density at the thick- target footpoints. 

Return-current energy losses can also produce a downward break in the X-ray spectrum 
(Sect. 5). The break energy depends on both the thermal structure of the plasma in the 
flare loop and on the non-thermal electron flux density distribution. These spectral modifi- 
cations and their evolution throughout flares provide an important test for the presence of 
initially un-neutralized electron beams and the return currents they must drive to neutralize 
them. Although RHESSI observations provide substantial information about the structure 
and evolution of flare spectra, only a lower limit on the electron flux density can usually 
be determined. Observations and analysis sufficiently accurate and comprehensive to ver- 
ify the presence of return current energy losses as the cause of a spectral break are yet to 
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be obtained. On the other hand, significant evidence exists (Sect. 5.3) indicating that return 
current losses do have an impact on flare hard X-ray emission. 

A thorough comparison of measured flare spectra with theoretical spectra computed from 
models incorporating collisional and return current energy losses (including their effect on 
the angular distribution of the non-thermal electrons), as well as nonuniform target ioniza- 
tion and albedo, is still needed. Spectral fitting alone, however, is not likely to distinguish the 
importance of these different mechanisms. Comparison of the time evolution of the spectra, 
as well as of the spatial structure of the X-ray emission, with expectations would certainly 
enhance the possibility of success for such an endeavor. 

The analysis of the evolution of X-ray source positions and sizes with photon energy and 
time provides another important test of the collisional thick-target model (Sect. 7). For these 
flares that show non-thermal source evolution in the corona and upper transition region, the 
source position and size are sensitive to the energy losses experienced by the non-thermal 
electrons. They are, in fact, sensitive to the very assumption that the sources are produced 
by electrons as they stream downward from an acceleration region higher in the corona. 
Further studies of the evolution of these coronal X-ray sources should substantially clarify 
the applicability of the collisional thick-target model. 

For completeness, we note that under some circumstances other radiation mechanisms 
may significantly contribute to the X-ray emission from non-thermal electrons. The pos- 
sibility that recombination (free-bound) radiation from the non-thermal electrons is some- 
times important is discussed in Brown and Mallik (2008, 2009) (also see Kontar et al. 2011). 
However, the contribution of non-thermal free-bound radiation has recently been found to 
be less significant than originally estimated (Brown et al. 2010). MacKinnon and Mallik 
(2010) have concluded that inverse Compton radiation may significantly contribute to the 
X-ray/y-ray emission from low-density coronal sources. 

Another testable aspect of the collisional thick-target model is the heating of the flare 
plasma by the non-thermal electrons. If the flare plasma is primarily heated by these elec- 
trons and the thick-target region is primarily in the chromosphere and lower transition re- 
gion, heating originating in the footpoints and expanding into the rest of the flare loop 
through “chromospheric evaporation” should be observed. On the other hand, if the loop is 
dense enough for the thick-target region to extend into the corona or if return-current heating 
is important, localized coronal heating and different ion abundances should be observed. 

It has generally been difficult to establish a clear connection between the location and 
evolution of X-ray sources produced by non-thermal electrons and by thermal plasma at 
different temperatures. This is largely because of a lack of high-cadence EUV images and 
spectra covering a broad range of coronal and transition region temperatures prior to the 
launch of SDO. Future studies of the coevolution of non-thermal X-ray sources and thermal 
sources in flares will be important in determining the extent to which heating mechanisms 
other than collisional heating by non-thermal electrons is significant. 

Predicting the expected evolution of the heated plasma is hampered by insufficient 
knowledge of the dominant heat transport mechanisms. We have seen evidence that many 
flares cool by classical thermal conduction or radiation once the heating has subsided 
(Sect. 8.4), but this is not likely to be the dominant transport mechanism during rapid heat- 
ing. Nevertheless, the spatial evolution of flare X-ray sources has so far been found to be 
consistent with chromospheric evaporation (Sect. 7.2). Also, the Neupert effect, observed 
in most flares, and Doppler-shift measurements qualitatively support the thick-target model 
(Sect. 8.3), but these do not rule out the possibility of other heating mechanisms temporally 
correlated with the electron beam collisional heating. As discussed in Sect. 3, substantial 
progress has been made in deducing the energy flux (total power) carried by non-thermal 
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electrons, but we usually can deduce only a lower limit to this energy flux. Continuing 
studies of flares similar to SOL2002-04- 15X03:55 (Ml. 2) and the initially cooler, early- 
impulsive flares (Sect. 3.5) may provide a better handle on this energy flux for comparison 
with thermal evolution. The thermal properties, energetics, and evolution of flares are dis- 
cussed further in Fletcher et al. (2011). 

12.2 Implications of X-ray Observations for Electron Acceleration Mechanisms and Flare 
Models 

In Sect. 10 we addressed the X-ray spectra of hard X-ray sources sometimes observed above 
the top of the hot loops or arcades of loops observed in flares. We reviewed results indicat- 
ing that the spectra are qualitatively, but not quantitatively consistent with expectations for 
electrons passing through a thin-target or quasi-thick-target region on their way to the thick- 
target footpoints of the flare loops. The apparent failure of these relatively simple models 
is probably a manifestation of the more complex above-the-looptop X-ray source structure 
revealed by RHESSI observations. 

Before RHESSI, time-of-flight delays in hard X-ray timing indicated that electrons were 
accelerated in a region somewhat above the looptops of the hot flare loops in most flares 
(Sect. 8.1). Also, cusps were observed at the top of flare loops by Yohkoh (e.g.. Sect. 10), 
indicating a magnetic connection to the region above the hot loops. 

RHESSI images have revealed flares with double coronal sources, one at or just above the 
top of the hot loops and the other at a higher altitude above the lower source. The centroid 
of the lower source is higher in altitude at higher X-ray energies, while the centroid of the 
upper source is lower in altitude at higher X-ray energies, indicating that energy release 
occurred between these coronal sources (Sui and Holman 2003; Sui et al. 2004; Liu et al. 
2008, 2009b). In one flare, the upper source accelerated outward to the speed of a subsequent 
coronal mass ejection. The white-light coronagraph on the Solar Maximum Mission (Webb 
et al. 2003), the Large Angle and Spectrometric Coronagraph (LASCO) and the Ultraviolet 
Coronagraph Spectrometer (UVCS) on SOHO (Ko et al. 2003; Lin et al. 2005), RHESSI (Sui 
et al. 2005b), and TRACE (Sui et al. 2006a) have all provided direct evidence for the presence 
of an extended, vertical current sheet above the hot flare loops and below the coronal mass 
ejection associated with eruptive flares. These and related observations are discussed further 
in Fletcher et al. (2011). 

These recent observations strongly support the “standard” model of eruptive solar flares, 
in which the hot flare loops build up below a vertical current sheet where inflowing mag- 
netic fields reconnect and a magnetic flux rope forms above the current sheet to become a 
coronal mass ejection (see Fletcher et al. 2011; Zharkova et al. 2011). The rate of electron 
acceleration has been observed to be correlated with the rate at which magnetic flux is swept 
up by the expanding footpoints of flare loops and with the rate of looptop expansion (Qiu 
et al. 2004; Sui et al. 2004; Holman 2005), indicating that the electron acceleration rate 
is correlated with the rate of magnetic reconnection. On the other hand, the observations 
also indicate that the rate of electron acceleration in the impulsive phase of flares is greatest 
before a large-scale current sheet or soft X-ray cusp is observed (Sui et al. 2008). 

Initially, when the electron acceleration rate is highest, the current sheet may be short and 
associated with slow-mode shock waves, as in Petschek reconnection. Fast reconnection jets 
(e.g., Wang et al. 2007) can stream upward and downward from the current sheet, possibly 
ending in fast-mode shock waves where they collide with slower magnetized plasma at 
the flare loop tops and the lower boundary of the magnetic flux rope (termination shocks). 
The pair of above-the-looptop X-ray sources may be associated with these fast-mode shock 
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waves. We have described possible evidence for these shock waves from radio observations 
in Sect. 11. 

The most difficult task is determining the dominant acceleration mechanism or mecha- 
nisms responsible for the energetic particles. The region above the flare loops contains or 
can contain quasi-DC electric fields, plasma turbulence, slow- and fast-mode shock waves, 
and collapsing magnetic traps, allowing for almost any acceleration mechanism imaginable. 
The problem is as much one of ruling out mechanisms as of finding mechanisms that work 
(cf. Miller et al. 1997). Acceleration mechanisms are addressed in Zharkova et al. (2011). 

In Sect. 9 we addressed the soft-hard-soft evolution of flare X-ray spectra. This spectral 
evolution could occur during the propagation of the electrons from the acceleration region 
to the thick-target footpoints. Return current losses, with their dependence on the electron 
beam flux (Sect. 5), for example, could be responsible for this evolution. However, the obser- 
vation that above-the-looptop sources also show this spectral evolution (Sect. 10.3) indicates 
that it is a property of the acceleration process rather than electron beam propagation. We 
saw in Sect. 9.2 that the soft-hard-soft behavior can be reproduced in the acceleration region 
if the acceleration or trapping efficiency first increases and then decreases. 

Flares displaying soft-hard-harder spectral evolution are of special interest, because they 
have been shown to be associated with high-energy proton events in space (Kiplinger 1995; 
Saldanha et al. 2008; Grayson et al. 2009). What is the connection between the appearance 
of energetic protons in space and X-ray spectral hardening late in flares? The answer to this 
question is important to both space weather prediction and understanding particle accelera- 
tion in flares. 

12.3 Implications of Current Results for Future Flare Studies in Hard X-rays 

What characteristics should a next-generation hard X-ray telescope have to make substantial 
progress in understanding electron propagation and acceleration in flares? The advances 
made with RHESSI have depended on its high-resolution count spectra that could generally 
be convolved with the detector response to obtain reliable photon flux spectra. These have 
been the first observations to allow detailed information about the evolution of accelerated 
electrons and associated hot flare plasma to be deduced for many flares. Equally important 
has been the ability to produce hard X-ray images in energy bands determined by the user 
during the data analysis process. This imaging capability has been critical to determining the 
origin of the X-ray emission at a given photon energy and in obtaining spectra for individual 
imaged source regions. These high-resolution imaging spectroscopy capabilities will remain 
important for continued progress. 

RHESSI’s X-ray imaging capability has allowed a clear spatial separation to be made 
for many flares between footpoint sources with non-thermal spectra at higher energies and 
looptop sources with thermal spectra at lower energies. However, in the energy range of 
overlap between ~10 keV and 50 keV, where both types of sources may coexist, it is of- 
ten difficult to distinguish weaker coronal sources (both thermal and non-thermal) in the 
presence of the stronger footpoint sources. This is because of the limited dynamic range of 
<100:1 (and significantly less for weaker events) that is possible in any one image made 
from RHESSI data. This is a consequence of the particular form of the Fourier-transform 
imaging technique that is used. Thus, in most flares the usually intense footpoints mask 
the much weaker coronal hard X-ray sources that can sometimes be seen in over-the- 
limb flares when the footpoints are occulted (e.g., Krucker and Lin 2008). In fact, these 
coronal hard X-ray sources can extend to high energies (up to ~800 keV, Krucker et al. 
2008b) and seem to be non-thermal in origin, thus making them of great interest in lo- 
cating and understanding the particle acceleration process. It is important to study these 
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non-thermal coronal sources in comparison with the footpoint sources, something that is 
currently not possible with RHESSI’s limited dynamic range except in the few cases with 
exceptionally strong coronal emission (see Sect. 10). In addition, again because of the 
RHESSI dynamic range and sensitivity limitations, it has not generally been possible to 
observe the thin-target bremsstrahlung emission from the corona that must be present from 
the electrons streaming down the legs of magnetic loops and also from electrons stream- 
ing out from the Sun and producing type-III bursts (see, however, Krucker et al. 2008c; 
Saint-Hilaire et al. 2009). For all of these reasons, a significantly greater dynamic range will 
be an important goal for future advanced solar hard X-ray instruments. 

Flares at the solar limb for which the hard X-ray footpoints are occulted by the disk pro- 
vide an important way of observing coronal hard X-ray sources, but these flares do not allow 
a comparison to be made between the coronal emission and the thick-target footpoint emis- 
sion. A possible substitute for a high-dynamic-range instrument is hard X-ray observations 
from two or more spacecraft. Under the right conditions, one spacecraft can observe all the 
flare emission while the other observes only the coronal emission, with the footpoint emis- 
sion occulted by the solar disk. Multi-spacecraft observations would also be important for 
deducing the directivity of the flare emission (especially in conjunction with X-ray polariza- 
tion measurements — see Kontar et al. 2011) and 3-D source structure. This multi-spacecraft 
approach, however, limits the number of flares for which the coronal and footpoint emissions 
can be compared. 

Hard X-ray timing studies have provided valuable information about electron propaga- 
tion and the location of the acceleration region (Sect. 8). Since the time of flight of energetic 
electrons from a coronal acceleration region to the thick-target loop footpoint is typically 
~10-100 milliseconds, the photon count rate must be high enough to distinguish differences 
in flux on these time scales. Time-of-flight studies have not been successful with RHESSI, 
because of its relatively low collecting area and, therefore, count rate. An instrument with 
the collecting area and pulse-pileup avoidance of CGI?0/BATSE, and the imaging and spec- 
tral resolution of RHESSI, would provide a new generation of studies on the characteristic 
time scales of propagation for the hard-X-ray-emitting electrons accelerated in flares. Alter- 
natively, smaller instruments sent closer to the Sun on, for example. Solar Orbiter or Solar 
Probe could achieve the required sensitivity. Flare studies on these time scales would provide 
important insights into the physical processes that impact the acceleration and propagation 
of energetic electrons in flares. 
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